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(N i ABSTRACT 

We present results of a sensitive 76 ksec Chandra observation of the young 
\q \ stellar cluster in NGC 2024, lying at a distance of ~415 pc in the Orion B giant 

molecular cloud. Previous infrared observations have shown that this remarkable 
cluster contains several hundred embedded young stars, most of which are still 
surrounded by circumstellar disks. Thus, it presents a rare opportunity to study 
\ X-ray activity in a large sample of optically invisible protostars and classical T 

Q-f Tauri stars (cTTS) undergoing accretion. Chandra detected 283 X-ray sources 

of which 248 were identified with counterparts at other wavelengths, mostly in 
the near-infrared. Astrometric registration of Chandra images against the Two 
Micron All Sky Survey (2MASS) resulted in positional offsets of ^0.25" near 
.^j | field center, yielding high confidence indentifications of infrared counterparts. 

The Chandra detections are characterized by hard heavily-absorbed spectra and 
spectacular variability. Spectral analysis of more than 100 of the brightest X-ray 
sources yields a mean extinction (A v ) ~ 10.5 mag and typical plasma energies 
(kT) ~ 3 keV. The range of variability includes rapid impulsive flares and persis- 
tent low-level fluctuations indicative of strong magnetic activity, as well as slow 
rises and falls in count rate whose origin is more obscure. Some slowly-evolving 
outbursts reached sustained temperatures of kT ~6 - 10 keV. Chandra detected 
all but one of a subsample of 27 cTTS identified from previous near and mid-IR 
photometry, and their X-ray and bolometric luminosities are correlated. We also 
report the X-ray detection of IRS 2b, which is thought to be a massive embed- 
ded late O or early B star that may be the ionizing source of NGC 2024. Seven 
millimeter-bright cores (FIR 1-7) in NGC 2024 that may be protostellar were 
not detected, with the possible exception of faint emission near the unusual core 
FIR-4. 
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Subject headings: open clusters and associations: individual (NGC 2024) — stars: 
formation — stars: pre-main-sequence - X-rays: stars 

1. Introduction 

The conspicuous HII region NGC 2024 lies about 15' east of the bright 09.51b supergiant 
( Orionis in the Orion B molecular cloud (= Lynds 1630). Its estimated distance is 415 pc 
(Anthony-Twarog 1982). Near-infrared (IR) observations have uncovered a remarkably dense 
cluster of young stars in the molecular cloud associated with NGC 2024 (Grasdalen 1974; 
Barnes et al. 1989 (B89); Comeron, Rieke, & Rieke 1996; Lada et al. 1991; Haisch, Lada 
& Lada 2000 (HLL); Haisch et al. 2001 (HLP); Meyer 1996 (M96)). The K-band survey of 
Lada et al. showed that the cluster covers an area of ^180 arc-min 2 and contains more than 
300 IR sources. The more sensitive survey of Comeron et al. detected ^150 sources in the 
central 30 arc-min 2 region of the cluster down to a completeness limit of K w 17, of which 
about two-thirds showed IR excesses characteristic of circumstellar disks. Using sensitive L- 
band observations, HLL found a somewhat higher infrared excess fraction >86%. A recent 
near-IR adaptive optics survey detected 73 stars, of which 3 are binaries and one is a triple 
system (Beck, Simon, & Close 2003 (B03)). Despite the high concentration of IR sources, 
there are no optically visible stars toward the cluster center. The cluster is undoubtedly 
young with age estimates ranging from ~0.3 My (Meyer 1996) up to a few My (Comeron et 
al. 1996). 

Millimeter and submillimeter observations provide evidence that star-formation is still 
underway in NGC 2024. Seven compact condensations of dust and gas denoted as FIR 1-7 
were detected at 350 - 1300 ^m by Mezger et al. (1988 (M88), 1992 (M92)), who argued that 
the condensations are isothermal protostars. These compact cores have been subsequently 
studied by numerous authors as summarized by Lai et al. (2002). Their evolutionary status 
is still a subject of debate, but the discovery of a faint near-infrared source at the FIR 
4 position (Moore & Chandler 1989; Moore & Yamashita 1995) and a unipolar redshifted 
outflow (Chandler & Carlstron 1996) suggests active star formation in the vicinity of FIR 
4. Polarization observations and Zeeman observations of absorption lines toward the cloud 
with radio interferometers have detected a magnetic field and there is some evidence that 
the field morphology around FIR 5 has been influenced by gravitational collapse (Crutcher 
et al. 1999; Lai et al. 2002). 

At X-ray wavelengths, ROSAT observations of NGC 2024 were analyzed by Freyberg & 
Schmitt (1995). The most sensitive image obtained was a 42 ks HRI exposure in which 52 
X-ray sources were detected. A shorter 22 ks PSPC image showed that most X-ray sources 
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have large absorption N H ~ 10 22 cm -2 and typical unabsorbed X-ray luminosities were L x 
~ 10 30 ergs s -1 . The ROSAT detections were well-correlated with known K-band sources 
and the X-ray properties suggested that most were embedded young stars. The ROSAT 
observations clearly revealed a large population of X-ray emitting sources in NGC 2024. 
But, because of the heavy X-ray absorption below ~1 keV and rather limited PSPC spectral 
resolution and bandpass («0.2 - 2.5 keV), spectral parameters for individual sources were 
not easily determined from ROSAT data. 

The NGC 2024 region is an ideal target for X-ray studies of star-formation for several 
reasons. It contains a diverse range of objects including class I protostars 1 , classical T 
Tauri stars (cTTS), weak-lined T Tauri stars (wTTS), the luminous infrared source IRS 2b 
which may be a massive young OB star, and seven millimeter condensations that may be 
protostellar. NGC 2024 lies at galactic coordinates (b,l) ~ (—16°, 207°) so source confusion 
from the galactic plane or galactic center is not an issue. In addition, the high absorption 
toward the cluster center reduces background contamination. 

We present here the results of a sensitive 76 ks Chandra observation of NGC 2024 
which provides higher angular and spectral resolution, broader energy coverage, and longer 
time monitoring than has previously been obtained. Our primary goals were: (i) to obtain 
a sensitive high angular resolution X-ray census of NGC 2024, (ii) to compare the X-ray 
population with existing IR catalogs and identify potential cloud members that might have 
escaped IR detection, (iii) to quantify the X-ray spectral and timing properties of the X-ray 
population, (iv) to determine if the total X-ray luminosity contributes significantly to the 
ionization of the visible HII region, and (v) to search for X-ray emission from the compact 
millimeter cores FIR 1-7. In the following we discuss the spatial, temporal, and spectral 
properties of the X-ray sources in NGC 2024 and compare the X-ray and infrared populations. 



2. Chandra Observation and Data Reduction 
2.1. Observation 

Chandra observed NGC 2024 from 8 August 2001 at 06:37 UT until 04:29 UT on 9 
August (Obsld 1878), spanning 21.9 hours. Exposures were obtained in faint data mode 
with a 3.2 s frame time using the ACIS-I imaging array as the primary detector. ACIS-I 



1 We adopt the the IR classification scheme used in Haisch et al. (2001) based on the IR spectral index a 
= d \og(XF\)/d logA evaluated in the 2.2 - 10 range. Class I sources have a > 0.3, flat spectrum sources 
have —0.3 < a < 0.3, class II sources have —1.6 < a < —0.3, and class III sources have a < —1.6. 
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consists of four front-illuminated 1024 x 1024 pixel CCDs with a pixel size of ~0.492" and 
a combined field-of-view (FoV) of «16.9' x 16.9'. The S2 and S3 CCDs in the ACIS-S array 
were also enabled, but our discussion here will focus on the ACIS-I data. The exposure 
livetime varied slightly between CCDs but was in the range 76,655 - 76,658 s. The Chandra 
aimpoint was at nominal RA = h h 41 m 46.3 s , DEC = -1° 55' 11.7" (J2000), which is ^32" 
from the strong millimeter source FIR-5 (M88). More detailed information on Chandra and 
its instrumentation can be found in the Chandra Proposer's Observatory Guide (POG) 2 . 

2.2. CIAO Data Reduction 

We followed standard data reduction procedures for ACIS data using CIAO software 
vers. 2.2. 1 3 . Our data reduction was based on the Level 1 files generated during standard 
processing by the Chandra X-ray Center (CXC) using vers. 2.6 of the Chandra calibration 
database. The CIAO data reduction included (i) application of the most recent aspect and 
focal length corrections, (ii) event selection based on good ASCA grades 0,2,3,4,6 and a clean 
status bit (status = 0) and removal of afterglow events, and (iii) energy filtering. For image 
analysis and source identification we used events in the 0.5 - 7 keV range to reduce particle 
background. Pixel randomization, which is introduced by CXC as part of the standard 
data processing to remove possible spatial aliasing effects, was retained. The algorithm to 
partially correct for the effects of charge transfer inefficiency (CTI) was not yet implemented 
in CIAO vers. 2.2.1 and we thus did not apply CTI corrections. 

2.3. Spatial Resolution and Limiting Sensitivity 

The on-axis spatial resolution of ACIS-I is limited by the physical pixel size of ~0.492". 
For on-axis point sources, approximately 90% of the encircled energy lies within 4 pixels 
(1.97") of the center pixel at 1.49 keV (POG). The limiting sensitivity of ACIS depends on 
many factors including the exposure time, background level, location of the source on the 
CCD-array, and intrinsic source properties. To obtain a rough sensitivity estimate for this 
observation, we consider an on-axis point source and adopt a 6 count detection threshold 
(or 7.83 x 10~ 2 c/ks for a livetime of 76.656 ks). We assume that the source has spectral 
properties typical of those sources detected in NGC 2024 (Sec. 4.3), namely an optically 



2 http:// cxc.harvard.edu/udocs/docs/docs.html 

3 Further information on Chandra Interactive Analysis of Observations (CIAO) software can be found at 
http://asc.harvard.edu/ciao/ 
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thin thermal plasma with kT ps 2.8 keV, absorption column density N H = 2.3 x 10 22 cm -2 , 
and abundances ~0.3 solar. For this typical source at an assumed distance of 415 pc, the 
unabsorbed luminosity detection limit from the Portable Interactive Multi-Mission Simulator 
(PIMMS) software is log L x ?a 28.8 ergs s -1 (0.5 - 7 keV). The detection limit is higher for 
off-axis sources, or for sources with heavier absorption or lower temperatures. Conversely, 
the detection limit is lower for sources having less absorption or higher temperatures. 

3. Data Analysis 

The X-ray data analysis consisted of (i) Chandra image generation and registration 
against IR positions from the 2MASS data base, (ii) identification of X-ray sources and 
their probable IR counterparts using existing catalogs, (iii) extraction and analysis of event 
lists for each source to determine source-specific quantities such as counts, mean photon 
energy, absorbed flux, and the Kolmogorov-Smirnov (KS) variability statistic, (iv) light 
curve extraction for all sources, and (v) extraction and fitting of spectra for brighter sources 
to obtain the absorption column density, characteristic plasma temperature, and unabsorbed 
luminosity. The analysis procedure is described in more detail below. 

3.1. Image Analysis and Source Identification 

We used the level 2 event file created by the CIAO processing described above to generate 
full- resolution 2800 x 2800 pixel images (0.492" pixels) as well as binned images. Light 
curves from on-chip background regions were inspected for large background fluctuations that 
might have resulted from solar flares and none were found. To insure accurate identification 
of infrared counterparts, we registered Chandra positions against IR positions in the the 
2MASS database 4 . To do this, we identified approximately 100 2MASS sources with obvious 
counterparts in the initial Chandra images. We then used the astrom program in the 
Starlink software package to derive a 4-coeflicient plate scale solution using 2MASS J2000 
positions and corresponding Chandra physical pixel positions. After applying the positional 
correction to the Chandra data, we found typical offsets between X-ray and 2MASS positions 
near the center of the ACIS-I detector to be a/ARA 2 + ADEC 2 « 0.25". These small offsets 
permitted IR counterparts to be identified with very high confidence. 

After applying the astrometric correction, we used the CIAO wavelet source detection 



http://www.ipac.caltech.edu/2mass/ 
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program wavdetect to search for X-ray sources. Several wavdetect passes were made 
using scale factors of 2,4, and 8 and false-alarm probabilities of 1CT 5 and 10~ 6 . The source 
lists from each wavdetect pass were compared and the Chandra images were visually 
inspected for missed or spurious detections. 

We identified 283 X-ray sources of which 248 were associated with known IR, optical, 
or radio counterparts based on searches of the SIMBAD data base, the 2MASS all-sky 
release Point Source Catalog, source lists from previous IR studies (HLL; HLPTL; B89; 
B03; M96), and a list of radio detections obtained in a recent 3.6 cm VLA survey conducted 
by Rodriguez et al. (2003). Table 1 lists the X-ray detections and their counterparts along 
with other source information. Table 1 is also supplied in electronic format. 

The positions and PSF-corrected 95% encircled energy regions around each source were 
used to extract source events. For the brightest sources a significant number of photons can 
fall outside the 95% encircled energy region, so larger 99% encircled energy regions were 
used for event extraction in bright sources. For on-axis sources the 95% and 99% encircled 
energy regions are nearly circular and have radii R95EE ~ 2" and R99EE ~ 8" (Feigelson et 
al. 2002). In the crowded cluster core region, it was necessary in some cases to use smaller 
extraction regions to avoid overlap from nearby sources. 

A background count rate of 0.17 counts s _1 per ACIS-I CCD in the 0.5 - 7 keV band 
(Table 6.6 of POG) was used along with the size of the source extraction region and exposure 
time to estimate the number of background counts for each source. The net source counts 
were then obtained by subtracting the number of background counts from the total number 
of counts in the source extraction region. The net counts for each source are given in Table 
1 along with other source information. The number of source counts that might have fallen 
outside the extraction region is within the net count uncertainty quoted in Table 1. Sources 
with fewer than 6 net counts were not included in Table 1. 

Figure 1 shows a smoothed broad-band Chandra image of a 4' x 4' region in the central 
part of the cluster where the X-ray source density is highest. Figure 2 is an unsmoothed 
image of the region in the vicinity of the infrared source IRS 2b, the possible ionizing source 
of NGC 2024 (Bik et al. 2003). Figure 3 is a zoomed smoothed broad-band image of the 
region near the millimeter sources FIR 1-7, which were not detected by Chandra with the 
possible exception of weak emission near FIR-4. This weak emission (4 counts) is shown in 
Figure 4. The astrometric-corrected positions of all 283 Chandra detections are plotted in 
Figure 5, overlaid with the positions of all 912 catalogued IR sources in the ACIS-I FoV. 
The X-ray source density is clearly higher near the center of the FoV, within the boundaries 
of the known IR cluster. 
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3.1.1. Unidentified Chandra Sources 

Figure 6 shows the positions of the 35 X-ray sources which lack counterparts at other 
wavelengths. About one-third of the unidentified sources have positions projected within 
the boundaries of the IR cluster and could be cluster members. However, the majority of 
the unidentified sources are probably extragalactic. 

The contribution of X-ray emitting foreground stars to the unidentified sample is ex- 
pected to be small because of the absence of optical counterparts (see also Barger et al. 
2002 for a breakdown of Chandra Deep Field North sources by object type). The faint hard 
emission of the unidentified X-ray sources also suggests a predominantly extragalactic origin. 
None of the 35 unidentified sources has more than 100 counts, while 43% of the identified 
sources have more than 100 counts. All but four of the unidentified sources have mean 
photon energies above the average computed for all X-ray detections (Sec. 4.2). 

Using the 2-8 keV X-ray number counts from Chandra Deep Field observations (Cowie 
et al. 2002), we estimate ~16 extragalactic sources in the ACIS-I FoV above our detection 
limit. This estimate should be treated with caution since (i) a reduction in the number of 
background sources detected toward NGC 2024 will occur as a result of absorption below a 
few keV from intervening molecular cloud material, which is non-uniform due to the north- 
south molecular ridge that transects NGC 2024, and (ii) the accuracy of the log N - log S 
distribution for extragalactic background sources from Chandra Deep Field observations at 
the lower Galactic latitude observed here (b m —16°) is not yet known, as already pointed 
out in a discussion of Chandra Orion Nebula observations (Feigelson et al. 2002). 

In Figure 6, we have circled six unidentified Chandra sources that could be extragalactic 
based on their exceptionally faint (<25 counts), hard ((E) > 4 keV), non- variable emission. 
Excluding these extragalactic candidates, there are seven remaining unidentified sources 
within the approximate IR cluster boundaries that may be young stars and for which follow- 
up observations at other wavelengths would be useful. These are sources 20, 21, 157, 168, 
180, 197, and 223. The latter source may have a 2MASS IR counterpart (Table 1). 

3.2. Timing Analysis 

To obtain a quantitative measurement of the variability for each source, we have com- 
puted the non-parameteric Kolmogorov-Smirnov (KS) statistic from the event list of each 
source (Table 1). The KS statistic is given by KS = y/n sup|/;(t) — fo(t)\ where n is the 
number of events, fi(t) is the normalized observed cumulative distribution and f (t) is the 
normalized model cumulative distribution, assuming a constant flux. A discussion of the KS 
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method and a comparison with the more familiar x 2 test is given by Rohatgi (1976) and 
additional information can be found in Babu & Feigelson (1996) and Press et al. (1992). 

If z is the value of the KS statistic then the null hypothesis probability P(z), or proba- 
bility of constant flux, is given by a convergent infinite series (Press et al. 1992): 



Truncating the series at six terms gives a few representative values: P(0.5) = 0.96, 
P(1.0) = 0.27, P(1.5) = 0.02, P(1.7) = 0.006, and P (2.0) = 6.7 x 10~ 4 . In the discussion 
below we adopt z = KS = 1.7 as the dividing line between variable and non- variable sources, 
where KS = 1.7 implies a variability probability P var = 0.994. That is, any source in Table 1 
with KS > 1.7 is considered to be variable. The KS values in Table 1 can be converted to the 
corresponding null hypothesis probability using the above equation. We note that the KS 
statistic is less reliable for diagnosing variability in weak sources (<25 counts) as discussed 
further below (Sec. 4.1). 

The KS statistic identifies those sources that were variable but does not provide infor- 
mation on variability profiles or timescales. To obtain this information, light curves were 
generated for all sources. Figure 7 shows light curves for ten sources in the ACIS-I FoV 
that demonstrate the wide range of variability detected. Figure 8 illustrates the continuous 
variability that was detected in source 94 throughout the observation. In Figure 9, we show 
the hard variable emission in source 207 with mean photon energies in excess of 4 keV and 
very high plasma temperatures. The interesting variability in these sources is discussed in 
more detail below (Sec. 4.5). 



We have used two different approaches to determine source spectral properties. First, 
the event list of each source was analyzed to determine the mean photon energy (E) and the 
absorbed flux F X)a b s (0.5 - 7 keV) values listed in Table 1. Second, spectra were extracted 
for the brightest sources and fitted with thermal plasma models to determine the absorption 
column density (N H ), plasma temperature (kT), and unabsorbed luminosity L x (0.5 - 7 keV) 
values given in Table 1. 



oo 




(1) 



3.3. Spectral Analysis 



These two methods are complementary. The advantage of event list analysis is that it 
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provides information on spectral hardness (or (E)) and F X)abs for all detections, including 
faint sources. Spectral fitting provides additional information on N H , kT, and L x but is only 
practical for brighter sources (>90 counts). 

The advantage of using both techniques is that the column density N H from spectral fits 
of brighter sources can be incorporated into the event list analysis. In that case, the event 
list analysis provides information on kT and L x , as well as (E) and F xa b s . If the event list of 
a specific source is first partitioned by time, then the time evolution of these four quantities 
can be plotted in "light curve" format, as shown in Figures 8 and 9. 

Event List Analysis: Our event list analysis followed the procedures described in Gagne, 
Daniel, & Skinner (2003), which we summarize here. The mean energy (E) is obtained from 
the event list, which contains the arrival time and energy of each photon. To compute the 
unabsorbed flux F xa b s , we generated an auxiliary response file (ARF) for each source, which 
contains information on effective area as a function of photon energy, taking into account 
the source position on the detector. The absorbed flux is then given by 

1 n E- 

where Ei is the energy of each photon, di(E) is the corresponding effective area, t is the 
livetime of the observation, A is the area of the source extraction region, and F bkg = 3.84 x 
1Q-18 er g g cm -2 s -i pixel- 1 i s the 0.5-7.0 keV background flux as determined from the CXC 
software tool make_acisbg. 

For brighter sources with >90 counts, additional information on the column density Nh 
is available from spectral fits (see below) and the event list analysis can be taken a step 
further to obtain plasma temperature T and L x . This is accomplished by generating a grid 
of simulated spectra using the IT vapec model in XSPEC vers. 11, with column densities 
in the range log N H = 20.0 - 24.0 cm~ 2 and temperatures in the range log T = 6.0 - 9.0 
(K) in increments of 0.1 dex. The grid is parameterized as a function of the mean photon 
energy (E), which is known from the event list of each source. Simulations show that (E) 
increases smoothly with Nh and T. The grid is then entered with the values of Nh and (E) to 
determine T. Once T is determined it can be used along with Nh to compute the conversion 
factor from absorbed to unabsorbed flux C(T,Nh) = F x /F Xja b s . The unabsorbed luminosity 
is then L x = 4vrd 2 F x = 47rd 2 F x , abs C(T,N H ). 

Spectral Fits: Spectra for sources with >90 counts were extracted along with source-specific 
response matrix files (RMF) and ARF files using the CIAO tool psextract. Background 
spectra were extracted from regions in the immediate vicinity of the source. Spectra of 
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weaker sources with ~90 - 300 counts were rebinned to 10 counts per bin prior to spectral 
fitting, while brighter sources (>300 counts) were rebinned to 15 counts per bin. We have 
restricted our spectral analysis to sources with >90 counts since our test fits show that the 
uncertainties in spectral parameters derived for weaker sources are often too large to permit 
definitive analysis. 

All sources with >90 counts were fitted with a single-temperature absorbed optically 
thin plasma model (IT vapec) using xspec vers. 11. A representative IT vapec fit of 
variable source 85 is shown in Figure 10, along with its pre-outburst spectrum. Brighter 
sources (>300 counts) were also fitted with a two-temperature 2T VAPEC model and in some 
cases the differential emission measure (DEM) model c6pvmkl was applied. During spectral 
fitting with vapec models, the free parameters were the equivalent neutral hydrogen column 
density (N H ), X-ray temperature (s) (kT), and normalization factors used to reproduce the 
observed count rate. 

Abundances were held fixed at solar values relative to Anders & Grevesse (1989), except 
for iron which was fixed at a lower value Fe = 0.3 x solar. The choice of a subsolar iron 
abundance was based on variable abundance fits of the brightest sources which invariably 
converged to a low Fe abundance near Fe ps 0.3 x solar. For example, we cite the relatively 
strong non- variable source 75 (3092 counts). Fits of this source with a IT vapec model 
converged to Fe = 0.22 [0.04 - 0.40] solar and 2T vapec models gave a very slight improve- 
ment with Fe = 0.35 [0.13 - 0.58] solar, where brackets enclose 90% confidence intervals. A 
comparison of IT vapec fits of different sources indicates that solar abundance models tend 
to overestimate Nh and unabsorbed L x relative to the improved fits obtained with subsolar 
Fe. 

The unabsorbed X-ray flux in the 0.5 - 7 keV range was measured from best-fit models 
after setting Nh = and then converted to an equivalent unabsorbed luminosity L x using a 
distance of 415 pc. Best-fit values of Nh and kT (with 90% confidence limits from the xspec 
steppar routine) and unabsorbed L x for the brightest sources are given in Table 1. 

We emphasize that kT values determined from the single-temperature IT vapec models 
are only used to define a characteristic time-averaged temperature that allows a large number 
of sources to be compared. This should not be construed to mean that the X-ray emitting 
plasma is isothermal. It is clear from our analysis that large temperature fluctuations did 
occur in some sources (Table 1 notes; see also Figs. 8 and 9). However, in sources with no 
statistically significant variability we only found definite improvement in spectral fits using a 
two-temperature model (2T vapec) in a few cases (Table 1 notes). This may be attributed 
to the strong absorption below ~1 keV that is present in most spectra, effectively masking 
any soft emission that could be present. The spectrum of source 85 in Figure 10 clearly 
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illustrates this low-energy absorption. 

3.4. Implications of Photon Pileup 

Pileup occurs when two or more photons are detected within the same pixel during one 
3.2 s readout period. Pileup is a consideration for high count-rate sources at small off-axis 
angles where the point-spread function (PSF) is narrowest, even for rapidly flaring sources 
where the count rate may increase to high levels only for short periods of time. In general, 
the effects of pileup are to make the spectrum appear harder than it actually is and to 
artificially reduce the count rate. However, there are other implications as discussed in the 
Chandra POG. 

We have identified those bright sources in NGC 2024 where significant pileup (>10%) 
could have occurred. Fortunately, pileup is not an issue for this observation except for a few 
sources, and then only at moderate levels. Three strongly variable sources located within k-V 
of the Chandra aimpoint had peak count rates that exceeded ~0.09 c s _1 and are suspected 
to have piled at the ~10% - 30% level for short periods. These sources are nos. 174, 192, and 
210. Another source which had a relatively high but steady count rate was no. 120, whose 
mean count rate was 0.10 c s _1 with short-term fluctuations up to ^0.15 c s -1 . However, it 
is offset ps2' from the aimpoint so its pileup should not have exceeded «10% - 15% due to 
off-axis broadening of the PSF. Because of moderate pileup in these four sources, their mean 
photon energies and flare temperatures given in Table 1 may be slightly overestimated, and 
total counts may be slightly underestimated. 

4. X-ray Properties 

We summarize here the global X-ray properties of the Chandra sources in NGC 2024 
based on temporal and spectral analysis. 

4.1. Time Variability 

A total of 72 out of 283 sources (25%) showed variability as inferred from KS > 1.7, and 
further information on sources with the strongest variability is given in the notes to Table 1. 
However, the variability fraction appears to be much larger among brighter sources. Only 
5% of those sources with less than 25 counts showed variability, whereas 64% of sources 
with more than 400 counts were variable. This trend is obvious in Figure 11, which shows a 



- 12 - 



monotonic increase in variability fraction with source counts. On physical grounds, there is 
little reason to suspect that the variability fraction should increase monotonically with the 
number of source counts. Instead, this trend can be attributed to the increased difficulty in 
detecting variability in weaker sources. It is thus likely that the true variability fraction of 
the 283 detections is greater than 25% due to undetected variability in faint sources, and 
this value should be interpreted as a lower limit. 

The range of variability shown in the light curves (Fig. 7) is truly amazing, consisting 
of rapid onset flares (e.g. sources 27, 81 and 174), variations with slow rises and decays (175, 
192, and 207), and sources with slowly declining or increasing count rates throughout the 
observation (94, 123, 269). Although rapid-onset flares are suggestive of solar-like behavior 
attributable to magnetic reconnection events, the origin of the slower variability is more 
obscure. In some cases the slow variability may be dynamically-induced, as for example 
by rotation of an X-ray bright active region on the star across the line-of-sight. Longer 
observations spanning days or weeks would be needed to search for periodic variability. 

4.2. Spectral Hardness (Mean Photon Energy) 

Figure 12 shows the distribution of mean photon energy (E) for all 283 detections. The 
average value is ( (E) ) = 2.7 keV. The distribution in Figure 12 is noticeably truncated at 
low energies below (E) ~ 1.4 keV. This is a result of the heavy absorption toward the cloud 
which effectively absorbs soft photons and skews the distribution. 

Figure 13 is a plot of (E) against the KS variability statistic. We have partitioned the 
diagram into four quadrants using KS = 1.7 to separate variable and non- variable sources 
and (E) = 2.7 keV to distinguish between sources with above average or below average mean 
photon energies. 

Several interesting features are apparent in Figure 13. First, all four quadrants are 
well-populated and there is no obvious relationship between variability (KS) and (E) for the 
sample as a whole. It is apparent that variability is not confined to hard sources since many 
sources with KS < 1.7 have (E) > 2.7 keV. The Kendall's r test in the ASURV statistical 
software package gives a probability of correlation between KS and (E) of p = 0.74 and the 
Cox test gives p = 0.93. Thus, statistical analysis does not provide conclusive evidence for 
a correlation. Second, it is obvious from Figure 13 that most of the 35 unidentified Chandra 
sources lie in the lower right part of the diagram and have hard non-variable emission. These 
sources are also faint (Sec. 3.1.1) and the hard emission suggests that most are extragalactic 
background sources but some could be deeply embedded sources in the cloud. Third, there 
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is a large spread in (E) among the Chandra sources, ranging from (E) = 0.92 keV (source 
181, which has the IR counterpart HLL-116) to the faint hard non- variable source 197 with 
(E) =4.77 keV and no known counterpart. 

4.3. Spectral Temperature and Absorption Estimates 

Since very few of the Chandra detections have optical counterparts, reliable information 
on the visual extinction Ay is generally not available. We have not made extensive use of 
methods for estimating Ay based on JHK photometry since these methods make a priori 
assumptions about dereddened colors of the star in order to extract the color excess and 
Ay. The errors which can be incurred in Ay are substantial, as discussed in more detail by 
Gagne et al. (2003). 

Figure 14 is a plot of Nh and kT derived from spectral fits of more than 100 of the 
brightest X-ray sources. It is apparent from Figure 14 that those sources with the highest 
absorption and highest temperatures are all variable. However, the converse is not true. 
Variability was also detected in sources that are situated in the lower left part of the diagram, 
where kT and Nh are relatively low. Thus, X-ray variability alone does not provide reliable 
information on either time-averaged temperature or absorption. 

The median absorption determined from 117 sources is log Nh = 22.18 cm -2 and the 
mean is log (Nh) = 22.37 cm~ 2 . The average temperature derived from IT vapec fits of 
sources with >90 counts is (kT) = 2.8 keV. In computing this average we have attempted 
to exclude time intervals during large rapid flares or slowly varying outbursts, although this 
was not possible in some cases where non-flare counts were low. 

The mean absorption column density log (Nh) = 22.37 cm -2 corresponds to a mean 
visual extinction (Ay) = 10.5 mag using the conversion Ay = 4.5 x 10~ 22 N H mag (Goren- 
stein 1975). This X-ray derived value is in excellent agreement with the value Ay ps 10.5 
mag determined for cluster members using IR data (HLL). We note a rather large spread in 
N H ranging from no detectable absorption in source 282 (= 2MASS 05422123-0159104) to a 
maximum of log N H = 23.18 cur 2 (Ay « 68 mag) in the variable source 160. 

Estimates of Ay based on IR data have been derived for 27 class II sources (cTTS) 
in NGC 2024 and the class III source (wTTS) HLP-24, as listed in Table 2 of HLP. We 
derived independent estimates of Nh for 21 of the class II sources from Chandra spectral fits 
(Table 1). Using the results of Gorenstein (1975) to convert from Nh to Ay, we generally 
find good agreement between the X-ray and IR determinations of Ay subject to rather large 
uncertainties in the X-ray estimates for fainter sources. But significant differences were found 
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for two sources. For source 31 (187 counts) whose counterpart is the IR source HLP-64, we 
obtain Ay = 3.1 [1.8 - 4.5; 90% conf.] mag while the IR data give Ay = 7.1 mag. For source 
69 (115 counts) the counterpart is HLP-14 which lies in a crowded IR field and the X-ray 
data yield Ay = 5.8 [3.6 - 8.6] mag but the IR data give a smaller value Ay = 0.8 mag. 

4.4. X-ray Luminosities 

Since reliable estimates of the visual extinction are not available for most of the em- 
bedded sources in NGC 2024, the intrinsic (unabsorbed) X-ray luminosities of many of the 
fainter Chandra detections are not yet known. We thus restrict our discussion of X-ray lu- 
minosities to the subsample of class II and flat-spectrum sources for which Ay estimates are 
available, and to an estimate of the total cluster X-ray luminosity. 

4-4- 1- X-ray Luminosities of Class II Sources (cTTS) 

NGC 2024 contains 27 class II sources identified from previous near and mid-IR ob- 
servations (Table 2 of HLP). Their spectral energy distributions peak in the near IR and 
decline more slowly toward the mid-IR than do normal stellar photospheres, indicative of 
circumstellar disks. 

Chandra detected 26 of the 27 class II sources listed in HLP. The only class II source not 
clearly detected was HLP-100 (= 2MASS 05413545-0152288). However, we do see a weak 
X-ray source (5 counts) at an offset of 0.75" from HLP-100, and we classify it as a possible 
detection. 

Figure 15 shows the cumulative X-ray luminosity distribution function for the 27 class 
II sources. This distribution function was obtained using the Kaplan-Meier (KM) estimator 
in the ASURV software package. We have estimated unabsorbed L x for the faintest class II 
detections using the PIMMS simulator and the absorbed fluxes in Table 1. The PIMMS 
simulations assumed a IT Raymond-Smith optically thin plasma model with a typical tem- 
perature for sources in NGC 2024 and an absorption N H computed from the Ay values given 
in HLP and the Gorenstein (1975) transformation. The possible detection of HLP-100 was 
treated as an upper limit (5 counts) with log L x (0.5-7 keV) < 28.63 ergs s _1 . 

A hard tail is clearly evident in the luminosity distribution function shown in Figure 
15. This is due to the anomalously high L x values of four class II sources with L x > 30.7 
ergs s -1 which exhibited rapid-onset flares or other prominent variability (source nos. 41, 
94, 133, 174). 
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Estimates of the bolometric luminosity (L bo i) were given for the 27 class II sources by 
HLP, ranging from L bo i = 61 L for HLP-74 (source 174) to L bo i = 1 L for both HLP- 
98 (source 220) and HLP- 100 (not detected). We used their L bo i values along with the 
unabsorbed L x (0.5 - 7 keV) measurements from Chandra to check for a correlation and to 
obtain a linear regression fit using the Buckley- James algorithm in ASURV. We excluded the 
four strongly variable sources noted above from the correlation test. The Kendall's r and 
Cox tests both give a correlation probability p = 0.98 between L bo i and unabsorbed L x (0.5 - 
7 keV). The regression analysis gives log L x (ergs s -1 ) = 29.86 + 0.54 (±0.19) log [L w /L ], 
with a standard deviation of 0.36 in log L x . For a typical class II source with L bo i ~ 1 L 
this relation yields log [L x /L 6o z] ~ —3.7. 

We consider the results of the above regression fit to be preliminary since it is based 
only on the IR-bright (K < 10.5) class II sources in NGC 2024 identified by HLP and lacks 
information on faint sources whose IR classification is not yet known. Even so, the regression 
fit is consistent at the ?slcr level with results from other young clusters having more complete 
optical/IR data such as IC 348 (Preibisch & Zinnecker 2002), and is nearly identical to the 
average (log[L x /L& z]) = —3.9 ± 0.7 found for a much larger sample of young stars in the 
Orion Nebula Cluster (Feigelson et al. 2002). Thus, based on this preliminary analysis, 
we find no reason to believe that the L x /L feoZ ratios of class II sources in NGC 2024 are 
substantially different from those in other young clusters. 

Previous studies of star-forming regions have reported correlations L x oc Lj oc L bo i, 
where Lj is the J-band luminosity (Casanova et al. 1995). Since the soft X-ray absorption 
in the energy range kT « 1.1 - 1.4 keV is comparable to that at J-band, it is possible to 
compare the observed (absorbed) soft-band X-ray and J-band luminosities. To compare with 
results for other star-forming regions, we have plotted the absorbed luminosity L xabs (1.1 
-1.4 keV) versus observed J magnitude from HLP for the known class II sources in NGC 
2024 (Fig. 16). Neither parameter has been corrected for extinction. We have excluded the 
four highly-variable class II sources noted above. There is a clear trend for brighter J-band 
sources to have larger L x a b s with the notable exception of one outlier, namely the bright 
IR source HLP-2 (Chandra no. 77). Excluding the outlier and the four strongly variable 
sources, the Cox and Kendall's r tests give a probability of correlation p = 0.99 between 
observed J magnitude and L xabs . Using the Buckley- James algorithm and again excluding 
the outlier, a regression fit gives log L Xjabs [l.l - 1-4 keV] (ergs s _1 ) = —0.37 (±0.1 1) (J— 12) 
+ 28.27 with a standard deviation of 0.56 in log L x ^ abs . Considering the uncertainties, the 
above regression slope —0.37 (±0.11) agrees quite well with the results for other star-forming 
regions (Casanova et al. 1995; Getman et al. 2002). 

It is noteworthy that several studies have shown obvious exceptions to the L Xjabs oc Lj 
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relation. In NGC 2024, the IR-bright source HLP-2 (J = 7.93, K = 6.56) is more than an 
order of magnitude less luminous in L xa b s than predicted by the above regression fit. We 
have undertaken an analysis similar to the above using K magnitudes and observed X-ray 
fluxes in the harder 3.2 - 4.2 keV band, for which the absorptions are comparable. Again, 
HLP-2 shows up as an outlier with a luminosity L X)a b s (3.2 - 4.2 keV) that is at least an order 
of magnitude below that expected for its Lk- Exceptions to the L Xja b s oc Lj relation were 
also seen in NGC 1333 (Fig. 11 of Getman et al. 2002) and the p Ophiuchi core (Fig. 4 
of Casanova et al. 1995). The exceptions that are bright at J-band but weak in X-rays 
such as HLP-2 and the objects SRI and SR3 in p Ophiuchi appear to have high bolometric 
luminosities and are thus more massive stars or perhaps unresolved multiple systems. 

X-ray Luminosities of Flat- Spectrum Sources 

Six flat-spectrum sources were identified in Table 1 of HLP, of which three were detected 
by Chandra. These were HLP-1 (source 187), HLP-12 (source 71), and HLP-38 (source 125). 
HLP-1 is the brightest K-band source in the HLP catalog (K fx 5.5 - 5.8). The IR-bright 
flat-spectrum source HLP-6 was not detected but it fell in a gap between CCDs. Although 
the sample is small, the 50% detection rate of these six flat-spectrum sources is lower than 
the 96% detection rate for known class lis. The X-ray luminosities of the three flat-spectrum 
detections were similar and in the range log L x = 30.14 - 30.37 ergs s -1 (Table 1). Because 
of the small sample size, statistical analysis of the flat-spectrum sources was not undertaken. 

4.4-3. Ionizing Source and Total Cluster X-ray Luminosity 

The mean absorbed X-ray flux of all 283 detections from Table 1 is (F xa b s (0.5 - 7 
keV)) = 4.6 x 10~ 14 ergs cm -2 s^ 1 . Assuming a IT optically thin plasma model with source 
properties typical of NGC 2024 (kT = 2.8 keV, N H = 2.3 x 10 22 cm -2 ), the corresponding 
mean unabsorbed flux is (F x (0.5 - 7 keV)) = 1.2 x 10~ 13 ergs cm -2 s _1 which equates to an 
unabsorbed luminosity (L x (0.5 - 7 keV)) = 2.5 x 10 30 ergs s" 1 . This average is about 1300 
times larger than the X-ray luminosity of the quiet Sun, underscoring the prodigious X-ray 
output of pre-main-sequence stars. Summing over 283 sources gives a total luminosity L x tot 
(0.5 - 7 keV) ~ 7 x 10 32 ergs s _1 . 

The source of the ionizing radiation which excites NGC 2024 has been difficult to identify 
but previous studies based on energy considerations have concluded that the exciting star is 
most likely of spectral type B0 - 09 (Hjellming 1968; Gordon 1969; Thronson et al. 1984). 
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Although the 09.51b star ( Ori may contribute to the ionization, it lies ~15' away and is not 
believed to be the primary excitation source (Gordon 1969; B89). The lack of an obvious 
ionizing source led Barnes et al. (1989) to suggest that the combined contribution from the 
embedded IR population, rather than a putative BO-09 star, might be the source of the 
ionization. It is thus relevant to ask if the total X-ray emission associated with the cloud 
population might be an important contribution to the ionization. Our Chandra observation 
shows that the known X-ray population could contribute L X)tot ~ 10 33 ergs s _1 toward the 
ionizing luminosity, assuming a distance of ~415 pc. The actual X-ray contribution could be 
higher if many deeply embedded sources within the cloud have escaped detection because of 
very high extinction. But, the above estimate of L xtot is much less than the Lyman continuum 
luminosity of a a BO - 09 star (Thompson 1984) and we conclude that the presently known 
X-ray population cannot by itself account for the required ionizing radiation. 

Recently, it has been suggested by Bik et al. (2003) that the infrared source IRS 2b is 
the long-sought ionizing source. As shown in Figure 2, this source lies ~5" west of B89-IRS 2 
at P. A. m —76° (Jiang, Perrier, & Lena 1984; Nisini et al. 1994). Bik et al. derive a spectral 
type of 08V - B2V for IRS 2b. IRS 2b was detected by Chandra (source 182) and its X-ray 
emission may be variable (KS = 1.41, P var = 0.96). Due to possible contamination from the 
nearby source 183 (see below), spectral fit results are somewhat uncertain. The absorption 
column density determined from IT vapec models gives A v = 14.9 [12.6 - 18.9; 90% conf.] 
mag using the Gorenstein (1975) conversion, but 2T models converge to an intrinsically 
softer spectrum with higher absorption (Table 1). 

Interestingly, the Chandra image reveals two nearby X-ray sources located ml. 7" - 2.9" 
north of IRS 2b (Fig. 2). These are the faint X-ray source 180 and the brighter source 
183, the latter being hard and variable. Interestingly, neither of these two X-ray sources 
have known IR counterparts, but source 183 has been detected at 3.6 cm in recent VLA 
observations (Rodriguez et al. 2003). The variability in source 183 suggests that it is an 
embedded young star and higher sensitive high-resolution searches for an IR counterpart are 
warranted. 



4.5. Interesting X-ray Sources 

Three Chandra sources showed exceptional variability that deserves further comment, 
namely sources 94, 192 and 207. 
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4-5.1. Persistent X-ray Variability in the cTTS Haro 5-59 

Source 94 is unusual in several respects. It was categorized as a class II source (cTTS) by 
HLP and is one of the few X-ray detections that has an optical counterpart, being associated 
with the emission-line star Haro 5-59 (V = 14.5). Its broad-band X-ray light curve shows a 
slowly declining count rate for most of the observation, with an abrupt increase near the end 
(Fig. 7). Closer inspection reveals that the variability is much more complex than suggested 
by Figure 7, as demonstrated in the four-panel light curves shown in Figure 8. 

Remarkably, the X-ray source is in a state of nearly continuous variability throughout 
the observation with significant variations (>3er) present in all four parameters: observed 
flux, (E), L x , and kT. The observed flux and L x generally undergo a slow decline which 
is punctuated by abrupt increases near UT 13 h and again about one hour prior to the 
end of the observation. The increases in flux and L x are accompanied by increases in mean 
photon energy (E) and kT. These fluctuations suggest reheating of the flaring plasma during 
the decay, but we cannot rule out the possibility that the fluctuations are due to flaring of 
different plasma structures - perhaps within the same active region. 

Figure 8 shows that the plasma temperature fluctuates from a low of kT ~ 2 keV 
up to kT ~ 5 keV, with the maximum value reached near the end of the observation. 
Because of these temperature variations, isothermal fits have difficulty in reproducing the 
time-averaged spectrum. More sophisticated fits of the time-averaged spectrum using a DEM 
model (c6pvmkl) show a double-peaked structure with maxima near ~1 keV and ~5 keV, 
confirming the multi-temperature nature of the plasma. 

The persistent X-ray variability in Haro 5-59 reflects strong surface magnetic activity 
and may be a consequence of nearly continuous low-level flaring that produces the energy 
input needed to maintain the plasma at temperatures above kT ~2 keV. This behavior is, at 
least qualitatively, similar to that recently detected in sensitive XMM-Newton observations 
of the nearby magnetically-active star Proxima Centauri (Giidel et al. 2003). The corona of 
Prox Cen displayed a nearly uninterrupted sequence of low-level X-ray flares, many of which 
were correlated with optical variability. Evidence for continuued reheating during flares has 
also been obtained in other active late- type stars such as Algol (Schmitt & Favata 1999). 
The variability detected here in Haro 5-59 may reflect similar processes occurring in a much 
younger magnetically active classical T Tauri star. 
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4-5.2. Slow Variability 

Prominent variability that lasted several hours was detected in sources 192 and 207. 
The variability is characterized by slow (non-impulsive) rise phases, as shown in Figure 7. 
These two sources exhibited the strongest X-ray variability of all detections as measured by 
their respective values KS = 27.4 and KS = 17.9 (Table 1). Source 192 is associated with 
the heavily reddened 2MASS source 05414611-0154147 (J = 17.30, H = 14.48, K = 11.42). 
Source 207 has a faint K = 14.8 IR counterpart (Meyer 1996). 

We have analyzed the time behavior of these two sources using light curves generated 
from time-binned event lists and by extracting time-partitioned spectra during the outbursts. 
In order to mitigate pileup, the central pixel in the PSF core was excluded in the spectra 
extracted for source 192 near the end of the observation when its count rate exceeded 0.1 c 
s _1 . The two analysis approaches give consistent results, as summarized below. 

The distinguishing feature of source 192 is its high temperature. Both analysis tech- 
niques indicate that the X-ray temperature reached very high values shortly after the be- 
ginning of the outburst, certainly in excess of kT ~ 6 keV and quite possibly as high as kT 
~ 10 keV. The temperature remained at values above kT ~ 4 - 5 keV during the four hour 
rise phase monitored by Chandra, and the observed X-ray flux continued to increase until 
the end of the observation. 

A similar situation exists for source 207, as shown in Figure 9. Remarkably, its temper- 
ature remained at levels above kT ~ 5 keV for most of the observation and may have risen 
to values in excess of kT ~ 10 keV during events such as that near UT 23 h. Spectral 
analysis confirms maximum temperatures of at least kT ~ 6 keV. 

In summary, sources 192 and 207 exhibit strong absorption (log Nh ~ 22.6 cm -2 , or Ay 
~ 18 mag) that is accompanied by spectacular variability and very high temperatures (kT > 
5-6 keV). Increases in (E) (hardening of the spectrum) and temperature increases appear 
to occur many hours after the beginning of the outburst. We have noticed similar behavior 
in Chandra observations of other young stars, such as the wTTS DoAr 21 in p Ophiuchus 
(Gagne et al. 2003). These results suggest that continued energy input may prolong the 
duration of some large outbursts. 

5. Comparison of X-ray and IR Properties 

We have used JHK magnitudes from 2MASS and other IR surveys to make compar- 
isons between the X-ray and IR properties of sources in NGC 2024. Figure 17 shows the 
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distribution of K magnitudes for those Chandra detections with known IR counterparts. The 
X-ray detections show a strong peak near K m 11 - 12. There is a strong tendency to detect 
the brightest K-band sources in X-rays and there is a rapid falloff in the number of X-ray 
detections for K > 12. Only 20% of the catalogued IR sources in the ACIS-I FoV with K > 
12 were detected in X-rays, while 74% of the IR sources with K < 12 were detected. 

The J-K versus K color magnitude diagram in Figure 18 shows that Chandra detections 
span a broad range of J-K colors, but no sources fainter than K = 16.7 were detected. The 
faintest IR source detected was source 108 in Meyer (1996) with K = 16.74, which corresponds 
to Chandra source 129. Figure 18 also reveals some notable non-detections. The class I source 
HLP-29 (K = 9.1) was not detected by Chandra, nor was the flat-spectrum source HLP-6 
(K = 7.6). But, the latter object lies in a gap between CCDs where sensitivity was reduced. 

The J-H versus H-K color-color diagram in Figure 19 allows us to discriminate between 
normally reddened stars and stars with true infrared excesses (as due, for example, to cir- 
cumstellar disks). The left and right dashed lines in Figure 19 show the normal reddening 
loci of MOV and AOV stars assuming intrinsic J-H and H-K colors from Bessell & Brett 
(1988) and the reddening law from Rieke & Lebofsky (1985). A significant number of the 
Chandra detections lie between the two dashed lines and are thus likely to be normally red- 
dened stars. More interesting are the sources lying to the right of the normal reddening band. 
These objects have true IR excesses and several lack IR classifications but were detected by 
Chandra. Thus, they are good candidates for new embedded cluster members. 

6. Candidate Protostars 

The NGC 2024 region contains the class I infrared source HLP-29 and seven compact 
millimeter sources FIR 1-7 which may be protostellar in nature (M88; M92). These objects 
represent early phases of stellar evolution about which little is known at X-ray wavelengths. 
We have thus carefully inspected the broad-band astrometric-corrected Chandra images for 
any evidence of X-ray emission from these candidate protostars. 

6.1. The Class I Source HLP-29 

No significant X-ray emission was detected from the class I source HLP-29 within a 
circular region of radius 3" centered on the J2000 position given by HLP. Only one photon 
was detected within the 95% encircled energy region centered at the IR position, which 
is consistent with the expected number of background events for a 76 ks exposure. Since 
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the extinction toward HLP-29 is not well-known, any conversion of the upper limit on its 
Chandra count rate to an equivalent upper limit on intrinsic L x would be highly uncertain. 

6.2. The Millimeter Sources FIR 1-7 

We searched for X-ray emission within 3" of the millimeter positions of FIR 1-7, including 
the (precessed) positions of FIR 1-6 from M88, FIR-7 from M92, and FIR 2-7 from the OVRO 
observations of Chandler & Carlstrom (1996). Because of the high absorption toward these 
dust condensations (M92), any X-ray photons that escape are likely to be hard. Thus, we 
used images over a slightly broader energy range of 0.5 - 8 keV (rather than 0.5 - 7 keV) 
to search for emission in the vicinity of FIR 1-7. Even so, no significant X-ray detections 
were found. However, we do find weak emission near FIR-4, as shown in Figure 4. Four 
photons were detected within 1.2" of the FIR-4 position given by M88, and three of these 
were localized in two adjacent pixels. The emission is hard and the mean energy of the four 
photons is (E) = 5.9 keV. Such emission would most likely be due to a flaring source deeply 
embedded within the cloud or an extragalactic source seen through the cloud. However, 
contamination from extragalactic sources along this line-of-sight is expected to be low due 
to the very high column density toward the FIR-4 dust condensation (M92). No more than 
one of the four photons is attributable to quiescent background in a 76 ks exposure (Table 
6.6 of POG). Furthermore, the probability of detecting two background photons in the same 
pixel is no larger than p m 2.6 x 10 -4 , assuming a maximum quiescent background rate in 
the 0.5 - 8 keV range of ^0.22 counts s _1 per ACIS-I CCD. The photon arrivals are well 
separated in time so the weak emission is not due to a cosmic ray afterglow event. The 
proximity of the hard X-ray emission to the FIR-4 millimeter peak and the unusual nature 
of this dust condensation warrant further comment. 

FIR-4 is exceptional in that it is the only one of the seven condensations that has a 
near-IR source at the position of the millimeter peak. A faint K-band source was discovered 
by Moore & Chandler (1989) and our inspection of 2MASS images confirms this K-band 
source at an offset of only 0.5" from the FIR-4 mm peak position of M88 (Fig. 4). The IR 
source was detected at H, K, and L bands by Moore & Yamashita (1995 = MY95). It is 
heavily reddened (H - K = 3.7) and there is a bright reflection nebula extending psIO" to 
its northwest. MY95 suggested that the reflection nebula is due to scattering in a region 
channeled out by a molecular outflow, and noted that the IR source might even be a localized 
peak in the surface brightness of scattered light rather than the actual illuminating source. 
The redshifted lobe of an outflow has been detected with the Owens Valley millimeter array 
by Chandler & Carlstrom (1996). 
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Even though the IR source is nearly coincident with the FIR-4 millimeter peak, the 
physical relationship between the two is not yet clear. It was argued by M92 that the 
IR source does not lie at the center of the FIR-4 dust condensation, since it would then 
be obscured by an extinction Ay ~ 1000 mag and be undetectable in the near-IR. The 
extinction of the IR source was estimated to be Ay < 50 mag by MY95. Similarly, any 
physical connection between the weak X-ray emission and either FIR-4 or the reddened IR 
source is not yet established since the X-ray peak is displaced sal" to the north (Fig. 4). 

7. Summary 

The embedded infrared cluster in NGC 2024 provides an interesting contrast to other 
well-studied star-forming regions in Orion such as the Orion Nebula Cluster (ONC). NGC 
2024 lies at approximately the same distance as the ONC but has a lower stellar density and 
is more heavily obscured, with a dearth of optically visible stars. Although star-formation 
appears to have already ceased in the ONC (Feigelson et al. 2002), active star- formation 
may still be underway in NGC 2024 and many of the cluster members are still in the class 
II accretion phase. 

The catalog of 283 Chandra sources is large enough for statistical studies and provides 
new information on the distribution of young stars within the cloud, absorption, X-ray 
spectral and variability properties, and energetics of the X-ray emitting population. The 
most important results to emerge from this study are the following: 

1. Chandra detected 283 sources in NGC2024 within a wl7' x 17' ACIS-I region centered 
near the millimeter source FIR-5, a five-fold increase over previous X-ray surveys. A 
total of 248 X-ray sources (88%) were identified with counterparts, mostly infrared 
sources. 

2. Roughly half of the 35 sources without counterparts are expected to be extragalactic 
based on source counts from Chandra Deep Field surveys. About a dozen of the 
unidentified sources have positions within the projected boundary of the known IR 
cluster and some of these could be embedded young stars. We have provided accurate 
X-ray positions based on registration of Chandra images with the 2MASS data base 
(Table 1) that will facilitate future searches for counterparts to these X-ray sources in 
the infrared and other spectral regions. 

3. The X-ray sources in NGC 2024 are in general characterized by hard heavily absorbed 
emission with strong variability. By fitting the spectra of more than 100 of the brightest 
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detections, we have derived a mean X-ray energy (kT) = 2.8 keV and mean extinction 
(Ay) = 10.5 mag. This latter value is in excellent agreement with extinction estimates 
from IR studies. However, considerable variation about these mean values was ob- 
served, with kT ~6 - 10 keV derived for some strongly- variable sources and extinctions 
as large as Ay ~ 68 mag for deeply embedded sources. 

4. Variability was detected at high confidence levels in 25% of the X-ray sources, but 
we have argued that the true variability fraction is likely to be larger because of the 
difficulty in detecting variability in weak sources. We find no compelling evidence 
that variability and spectral hardness are correlated. Of particular interest are sources 
which undergo poorly-understood slow variability lasting many hours during which 
plasma temperatures are maintained at high values kT > 5 - 6 keV such as sources 
192 and 207 (Fig. 7). Such behavior is decidedly different from solar-like impulsive 
bursts (e.g. sources 27, 81 and 174) that signal magnetic reconnection and rapid 
energy release. Our time-resolved spectroscopic analysis suggests that the plasma in 
such slowly varying sources is being reheated. However, we cannot spatially resolve 
individual flaring regions so it is possible that multiple unresolved flares are taking 
place, perhaps in the same active region. We also find clear evidence for persistent 
low-level variability in the star Haro 5-59, which is probably the consequence of nearly 
continuous low- level flaring in this magnetically active cTTS. 

5. A comparison with IR data shows that the brightest K-band sources are preferentially 
detected in X-rays. Chandra detected 74% of the known IR sources with K < 12 but 
only 20% of those with K > 12. Near-infrared color-color diagrams reveal that many 
Chandra detections are heavily reddened, including extreme objects such as the hard 
variable X-ray source 192 (= 2MASS J0541461 1-0154147; J-H = 2.82, H-K = 3.06; Fig. 
19). These heavily reddened sources could be young cluster members and follow-up 
observations in the mid-IR would be useful to determine if disks are still present. 

6. The detection rate of 27 known class II sources (cTTS) in NGC 2024 is at least 96%, 
and may be as high as 100%. The X-ray and bolometric luminosities of these class II 
sources are correlated and we also find a correlation between the observed luminosities 
in the equally absorbed J band and 1.1 - 1.4 keV X-ray band. A notable exception 
is the IR-bright cTTS HLP-2, whose X-ray emission is at least an order of magnitude 
lower than expected from regression fits. 

7. The ionizing source of NGC 2024 has been elusive, but previous studies have determined 
that the required ionization flux is equivalent to that of a BO-09 star. We estimate the 
total X-ray luminosity from all 283 X-ray sources in the Chandra FoV to be L X)to t (0.5 
- 7 keV) ~ 10 33 ergs s" 1 at an assumed distance of ~415 pc. This is significantly less 
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than the Lyman continuum luminosity of a BO-09 star. Thus, the X-ray contribution 
from the known population cannot account for all of the missing ionizing flux. Recent 
work by Bik et al. (2003) suggests that the luminous IR source IRS 2b is a young 08V 
- B2V star and may be the long-sought ionizing source. IRS 2b is detected by Chandra 
and its X-ray emission is likely variable, possibly a sign of magnetic activity in this 
embedded high-mass star. 

8. Close examination of the Chandra images shows no significant X-ray detections (>6 
counts) of any of the seven millimeter-bright condensations FIR 1-7 in NGC 2024. 
These condensations, which may be protostellar, are thus not yet sufficiently evolved 
to produce X-rays or any intrinsic X-ray emission has been completely absorbed by 
intervening material. A potentially interesting exception is FIR-4, for which we detect 
faint hard localized emission (4 counts) within 1.2" of the mm position. Previous 
observations have detected a heavily-reddened near-IR source within 0.5" of FIR-4 
as well as a CO outflow, indicating that stars have recently formed and may still be 
forming in this region. However, the physical relationship between the FIR-4 dust 
condensation, near-IR source, and faint X-ray emission is not yet clear because of the 
small positional offsets between the millimeter, IR, and X-ray peaks (Fig. 4). 
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Fig. I. — Smoothed broad-band (0.5 - 7 keV) Chandra ACIS-I image of the central w4' x 
4' region of NGC 2024 where X-ray source density is highest. Smoothing was done with the 
wavelet smoothing algorithm in the HEASARC XIMAGE software package (Rosati, Burg, 
& Giacconi 1994). Source numbers corresponding to the first column of Table 1 are shown 
for a few of the brightest sources. Crosses mark the peak positions of the compact millimeter 
cores FIR 1-6 from M88 (FIR 1-6) and M92 (FIR 7). Coordinates are J2000. 

Fig. 2. — Unsmoothed broad-band Chandra ACIS-I image of the region near IRS 2b (source 
182), which has been proposed as the ionizing source of NGC 2024 (Bik et al. 2003). Source 
183 lies 1.7" north of IRS 2b and has no known IR counterpart but has been detected with 
the VLA (Rodriguez et al. 2003). Source 180 has no known counterpart. Source 187 (B89- 
IRS 2) is a bright near-IR source with K = 4.6 - 5.8. The pixel size is 0.492" and coordinates 
are J2000. 

Fig. 3. — Smoothed broad-band Chandra ACIS-I image of the region in the vicinity of the 
millimeter cores FIR 1-7. The smoothing algorithm is the same as in Figure 1. Crosses mark 
the positions of cores FIR 1-6 (M88) and FIR 7 (M92). Open circles mark the positions of 
FIR 2-7 from the OVRO observations of Chandler & Carlstrom (1996). None of the cores 
were detected in X-rays with the possible exception of weak emission naear FIR-4 (see Fig. 
4). Two Chandra sources (nos. 122 and 139) in the vicinity of FIR 1 and FIR 2 are identified. 
Coordinates are J2000. 

Fig. 4. — Unsmoothed Chandra ACIS-I image (0.5 - 8 keV) showing 4 X-ray photons within 
1.2" of the mm source FIR 4. Pixel size is 0.492". Dark pixels contain two photons and light 
pixels contain one photon. The cross (+) marks the position of the mm/sub-mm source 
FIR-4 from the IRTF observations of M88, and a dashed circle of radius 1" centered on 
FIR-4 is shown for scale. The small open circle marks the position of the 3 mm continuum 
peak of FIR-4 determined with the Owens Valley millimeter array (Chandler & Carlstrom 
1996). An asterisk is at the position of the K-band source 2MASS 05414407-0154448 first 
identified by Moore & Chandler (1989). The three X-ray photons at the top of the dashed 
circle have energies 2.49, 7.39, and 7.77 keV while the photon just below FIR 4 has energy 
5.98 keV. Coordinates are J2000. 

Fig. 5. — Positions of all 283 Chandra detections in the ACIS-I field of view (filled circles), 
overlaid with positions of known IR sources (open circles). The IR source positions were 
taken from the all-sky release of the 2MASS Point Source Catalog and other surveys (B89, 
B03, HLL, HLP, M96). Coordinates are J2000. 

Fig. 6. — Positions of the 35 Chandra detections in the ACIS-I field of view without known 
counterparts at IR, optical, or radio wavelengths. Six circled sources have faint hard non- 
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variable X-ray emission and may be extragalactic (AGNs). The dotted line marks the ap- 
proximate boundary of the embedded IR cluster. Coordinates are J2000. 

Fig. 7. — Broad-band Chandra light curves of ten sources showing interesting variability. 
Source numbers correspond to the first column of Table 1. The count rate of source 192 
may be slightly underestimated toward the end of the observation due to moderate pileup 
at count rates above 0.1 c s _1 . 

Fig. 8. — Chandra light curves of variable source no. 94 (= J054138. 9-015936), whose optical 
counterpart is the emission-line star Haro 5-59. The panels show, from top to bottom, the 
absorbed flux f x (0.5 - 7 keV), mean photon energy (E), unabsorbed luminosity L x (0.5 - 
7 keV), and plasma temperature kT. Error bars are la. Note that the source is in a state 
of continuous variability during the ~21 hour observation, and the plasma may have been 
reheated at ^13 h and again at ^27 h. 

Fig. 9. — Chandra light curves of the hard variable source no. 207 (= J054147.4-015526) 
showing the same quantities as in Fig. 8. Error bars are la. 

Fig. 10. — Chandra ACIS-I spectra of source 85. Top panel: Data points with la error 
bars show the background-subtracted flare spectrum, binned to a minimum of 15 counts per 
bin. Dotted line is the background-subtracted pre-flare spectrum binned to a minimum of 
10 counts per bin (error bars omitted for clarity). Solid line is a best-fit IT vapec model of 
the flare spectrum using solar abundances except for iron, which was held fixed at Fe = 0.3 
solar. The best-fit model of the flare spectrum converges to an absorption N H = 1.8 [1.7 - 
2.0] x 10 22 cm" 2 and kT = 3.8 [3.3 - 4.3] keV with x 2 /dof = 85.8/104, where brackets enclose 
90% confidence intervals. The pre-flare temperature is kT 3.0 keV and the observed flux 
increased by about a factor of two during the flare. Bottom panel: Residuals of the best-fit 
IT vapec flare model, in the sense of data — model. 

Fig. 11. — Histogram showing the fraction of variable sources (KS > 1.7) as a function of 
the number of source counts. The apparent monotonic increase in variability fraction with 
source counts is very likely a selection effect caused by the inherent difficulties in detecting 
variability in weak sources. 

Fig. 12. — Histogram showing number of X-ray sources as a function of mean photon energy 
(E). The distribution is truncated at low energies below ssl.4 keV due to the absorption of 
soft photons by intervening cloud material. The mean for all detections is (E) = 2.7 keV. 

Fig. 13. — Distribution of mean photon energy (E) versus the KS variability statistic for 
all 283 X-ray detections in NGC 2024. Variable sources have KS > 1.7 and lie above the 
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solid line. The vertical dashed line is the average value (E) = 2.7 keV for the entire sample. 
Most Chandra detections without counterparts (open squares) are non-variable and have 
above-average (E). All four quadrants are well-populated and there is no obvious correlation 
between KS and (E). 

Fig. 14. — Absorption column density N H versus time-averaged plasma temperature kT as 
determined from IT vapec fits of the spectra of the brightest Chandra detections. For 
clarity, error bars (90% confidence) are shown only for non-variable sources. Circled sources 
showed strong variability in their light curves. Sources with the highest plasma temperatures 
and largest absorptions are all variable. But, the converse is not true since many sources 
with low temperature and absorption also varied. Thus, variability alone is not a reliable 
indicator of temperature or absorption. 

Fig. 15. — Kaplan-Meier estimator showing the cumulative X-ray luminosity distribution for 
known class II sources (classical T Tauri stars) in NGC 2024. The hard tail above log L x ps 
30.7 ergs s -1 is due to the contribution of four class II sources that experienced rapid-onset 
flares or other strong variability (sources 41, 94, 133, 174). 

Fig. 16. — Absorbed X-ray luminosity (1.1 - 1.4 keV) versus observed J magnitude for known 
class II sources in NGC 2024. No extinction corrections have been applied. In this energy 
range the X-ray and J-band extinctions are approximately equal. The source shown as an 
X-ray upper limit is HLP-100, which was possibly detected (see text). Four class II sources 
that experienced flares or other strong variability are omitted in this figure (sources 41, 94, 
133, 174). 

Fig. 17. — Histogram showing the number of Chandra detections versus K magnitude for 
those X-ray sources with known IR counterparts. The dotted line shows the corresponding 
histogram for all 912 known IR sources in the Chandra field-of-view. 

Fig. 18. — Color-magnitude diagram for IR sources within the Chandra ACIS-I field-of-view 
having J and K photometry. IR sources detected by Chandra are shown with filled symbols 
and undetected IR sources with open symbols. Undetected IR sources in the field which 
lack classifications based on IR spectral index are not shown. The class I source HLP-29 
(circled cross) was not detected. The undetected flat-spectrum source HLP-6 (K = 7.6) lies 
in a gap between CCDs where sensitivity is reduced. The faintest K-band source detected 
by Chandra is source 129, whose IR counterpart is source 108 of Meyer 1996 (M-108; K = 
16.74). Chandra source no. 160 (= M-165) is not visible in this plot and has K = 12.31, J-K 
= 8.59. 

Fig. 19. — Color-color diagram for IR sources within the ACIS-I field-of-view detected by 
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Chandra. Also shown are the colors of the undetected class I source HLP-29 (circled cross) 
and the undetected class II/flat-spectrum sources listed in HLP. Chandra detections are solid 
symbols and undetected IR sources are open symbols. Undetected IR sources in the field 
which lack classifications based on spectral index are not shown. The dotted line at lower 
left is the main sequence. The sloping dashed lines mark the approximate reddening band 
for normally reddened MOV (left) and AOV (right) stars, based on data from Bessel & Brett 
(1988) and Rieke & Lebofsky (1985). The slopes of the dashed lines match that of the A v 
= 10.5 mag reddening vector, which corresponds to the mean extinction toward NGC 2024. 
Sources to the right of the band have excess reddening. The most heavily reddened Chandra 
detections are identified. 
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Table 1. Chandra X-ray Sources in NGC 2024 a 



No. 


Name 


R.A. 


Decl. 


Net 


KS 


E 


N H 


kT 


F x ,abs 


log L x 


Identification 






(J2000) 


(J2000) 


Counts 




(kcV) 


(10 22 cm" 2 ) 


(keV) 


(erg/cm 2 /s) 


(erg/s) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



1 


J054112. 9-015459 


5 


41 


12.97 


-1 


55 


00.0 


29 ± 7 


1.25 


3.03 






5.33e-15 






2 


J0541 13.2-015332 


5 


41 


13.29 


-1 


53 


32.0 


57 ± 9 


0.80 


1.29 






5.33e-15 




05411324-0153306 


3 


J054114.4-015121 


5 


41 


14.43 


-1 


51 


21.0 


25 ± 7 


1.00 


2.14 






2.97c-15 




05411449-0151204 


4 


J054118. 0-014929 


5 


41 


18.00 


-1 


49 


29.7 


41 ± 8 


1.21 


2.74 






7.79c-15 




05411799-0149292 


5* 


J054118. 3-014953 


5 


41 


18.36 


-1 


49 


53.2 


120 ± 12 


0.88 


2.88 


1 O+0.8 




2.45e-14 


30.21 


05411838-0149531 


6 


J054119.1-015226 


5 


41 


19.14 


-1 


52 


26.9 


40 ± 8 


0.93 


2.00 






5.09c-15 




05411918-0152268 


7 


J054120.6-015857 


5 


41 


20.68 


-1 


58 


58.0 


545 ± 25 


2.08 


1.42 


2+° 1 


i 9 +0.1 
1- -0.3 


4.92c-14 


30.16 


Haro 5-46, 05412068-0158582 


8 


J054121.4-015629 


5 


41 


21.42 


-1 


56 


29.0 


17 ± 6 


1.08 


3.97 






7.18c-15 






9 


J054122.5-015041 


5 


41 


22.50 


-1 


50 


41.8 


41 ± 8 


0.53 


3.50 






1.20e-14 






10 


J054124.2-014934 


5 


41 


24.23 


-1 


49 


34.7 


666 ± 27 


0.68 


1.16 


04 +o os 


n S+ - 1 


5.47c-14 


30.09 


HLL-2,05412424-0149348 ^ 


11 


J054124.8-015203 


5 


41 


24.83 


-1 


52 


03.9 


8 ± 5 


0.81 


3.74 


3.76c-15 




05412475-0152063 ^ 


12 


J054125.0-015226 


5 


41 


25.01 


-1 


52 


26.9 


76 ± 10 


0.70 


2.17 






1.08e-14 




HLL-31,05412501-0152267 ^ 


13 


J054125.4-014809 


5 


41 


25.43 


-1 


48 


09.3 


46 ± 8 


0.60 


2.34 






5.68c-15 




05412553-0148092 


14 


J054125.5-020111 


5 


41 


25.55 


-2 


01 


11.1 


59 ± 9 


0.90 


1.31 






3.92c-15 




05412562-0201116 


15* 


J054125. 8-015727 


5 


41 


25.87 


-1 


57 


27.7 


25 ± 6 


0.69 


2.09 






4.00c-15 
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Uo414ooo-Uloo2zU,M-loo,Boy-ol 


155 


Tn tr ,1 1 ,10 o n ~i crooo 

JU5414d.8-U15.M8 


5 41 


43.85 


-1 


53 


38.3 


19 it 5 


0.88 


3.03 






3.43e-15 




Uo4143o4-U15oooo,M-lo4 


156 


Tntr-ii /io n n "i croon 

JU5414d.9-U15.329 


5 41 


43.91 


-1 


53 


29.1 


187 it 15 


1.42 


3.29 


9 v +0.8 


A n + 6.0 

4 ' 2 -1.5 


4.16e-14 


30.39 


0541439 l-U15ozyl, M-1 oo 


157 


J054143. 9-015554 


5 41 


43.94 


-1 


55 


54.8 


10 ± 4 


0.87 


3.65 






2.17e-15 






158 


J054143. 9-015845 


5 41 


43.96 


-1 


58 


45.1 


22 it 6 


2.01 


3.04 






4.12e-15 




HLL-139, 05414398-0158454 


159 


T0~>41 43 9-01 ~>302 


5 41 


43.99 


_l 


53 


02.2 


9 it 4 


0.62 


3.25 






1 95e-1 






160*f 


J054144.1-015347 


5 41 


44.12 


-1 


53 


47.1 


301 it 18 


6.24 


4.16 


i5.it«:J 


2 0+ 13 


8.54e-14 


31.38 


05414412-0153472,M-165 


161 


J054144.1-015520 


5 41 


44.17 


-1 


55 


20.4 


34 ± 7 


0.85 


3.82 






9.45e-15 




M-167 


162 


J054144. 1-015735 


5 41 


44.18 


-1 


57 


35.5 


8 it 4 


0.54 


1.87 






6.53e-16 




05414420-0157353 


163 


J054144.2-015404 


5 41 


44.22 


-1 


54 


04.9 


48 it 8 


1.07 


3.54 






1.14e-14 




05414422-0154051.M-170 


164 


J054144.3-015524 


5 41 


44.33 


-1 


55 


24.8 


34 ± 7 


0.58 


2.27 






4.21e-15 




M-173 


165 


J054144.3-015420 


5 41 


44.36 


-1 


54 


20.2 


74 ± 10 


1.48 


2.95 






1.34c- 14 




05414436-0154203,B-84,M-174 


166 


J054144.3-015335 


5 41 


44.37 


-1 


53 


35.8 


97 ± 11 


1.32 


3.44 


5.0l 2 ;? 


9 o+3.4 


2.19e-14 


30.30 


05414436-0153358,M-176 


167 


J054144.4-015522A 


5 41 


44.40 


-1 


55 


22.8 


138 it 13 


0.85 


2.44 


9 1+0.5 


9 c+0.9 


2. 07c- 14 


30.11 


HLL-107,HLP-58,05414440-0155229,B89-3 


168 


J054144.4-015522B 


5 41 


44.49 


-1 


55 


22.0 


56 it 9 


0.60 


3.28 




1.28e-14 






169 


J054144.5-015740 


5 41 


44.50 


-1 


57 


40.2 


266 it 17 


1.13 


1.99 


o-8±g;l 


2 1+ ' 4 


3.09e-14 


30.09 


HLL-133,05414451-0157402,B89-29 


170 


J054144.7-015348 


5 41 


44.74 


-1 


53 


48.5 


39 ± 7 


1.83 


3.91 






1.06e-14 




M-182,05414474-0153485 


171 


J054144.8-015425 


5 41 


44.82 


-1 


54 


25.1 


47 it 8 


1.04 


3.11 






8.77c-15 




HLL-77,HLP-22,05414482-0154251,B-85,V- 


172 


J054144.8-015435 


5 41 


44.84 


-1 


54 


35.9 


15 it 5 


0.79 


3.29 






3.08e-15 




05414483-0154357,B-86,M-184,HLP-73 


173 


J054144.9-015717 


5 41 


44.96 


-1 


57 


17.3 


41 ± 7 


0.55 


0.96 






2.59c-15 




05414496-0157174 


174*f 


J054145.0-015406 


5 41 


45.06 


-1 


54 


06.2 


651 it 27 


16.0 


3.72 




6 1+ 3 1 


1.70c- 13 


31.01 


HLP-74,05414506-0154063,M-188 
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No. 


Name 


R.A. 


Dccl. 


Net 


KS 


E 


N H 


kT 


Fx.abs log L x 


Identification 






(J2000) 


(J2000) 


Counts 




(keV) 


(10 22 cm" 2 ) 


(keV) 


(erg/cm 2 /s) (erg/s) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) (11) 


(12) 



175*f 


J054145.0-015144 


5 


41 


45.08 


-1 


51 


44.3 


1810 ± 44 


10.1 


3.03 


2 


o+0.2 

— 0.2 


3 


7+0.7 

— 0.5 


3.70e-13 


31 


26 


HLL-19,05414508-0151443 


176 


J054145. 2-015422 


5 


41 


45.23 


-1 


54 


22.7 


196 ± 15 


1.40 


3.03 


3 


9+1-2 

— 0.9 


1 


q+0.8 

— 0.5 


3.75c-14 


30 


59 


B-89,M-190 


177 


J054145. 3-015502 


5 


41 


45.33 


-1 


55 


02.6 


16 ± 5 


2.11 


3.00 










5.90e-15 






05414532-0155025, M-193 


178 


J054145. 3-015156 


5 


41 


45.38 


-1 


51 


56.6 


224 ± 16 


0.49 


3.22 


3 


4+1.1 

^-0.8 


2 


n+1.8 

y -0.9 


4.98e-14 


30 


52 


05414538-0151566 


179 


J054145. 3-015355 


5 


41 


45.40 


-1 


53 


55.7 


13 ± 5 


0.48 


2.58 








2.17c-15 






M-195 


180 


J054145. 4-015425 


5 


41 


45.42 


-1 


51 


25.7 


64 ± 9 


0.76 


3.97 










1.74c-14 








181 


J054145. 4-015616 


5 


41 


45.48 


-1 


56 


17.0 


8 ± 4 


0.19 


0.92 










3.96c-16 






HLL-116,05414548-0156169,M-196,B89-22 


182* 


J054145.4-015428 


5 


41 


45.48 


-1 


54 


28.6 


371 ± 20 


1.41 


3.14 


3 


O+0.9 
°-0.5 


2 


g+0.7 
— 0.6 


7.42c-14 


30 


65 


05414550-0154286,B-90,M-197,IRS2b,V-15 


183*f 


J054145.5-015426 


5 


41 


45.52 


-1 


54 


26.9 


220 ± 16 


5.94 


3.41 


2 


n+1.0 
9 -1.0 


6 


7+V-, 

— 6.1 


5.28e-14 


30 


39 


V-16 


184 


J054145.6-015732 


5 


41 


45.61 


-1 


57 


33.0 


34 ± 7 


0.95 


1.60 










3.47e-15 






05414562-0157330 


185 


J054145. 6-015504 


5 


41 


45.69 


-1 


55 


04.6 


14 ± 5 


0.75 


3.98 










6.56e-15 






M-200 


186 


J054145. 7-015212 


5 


41 


45.72 


-1 


52 


12.8 


10 ± 4 


0.57 


1.91 










8.49c-16 






05414576-0152130 


187* 


J054145. 7-015429 


5 


41 


45.80 


-1 


54 


29.8 


162 ± 14 


1.30 


3.42 


3 


i+l.O 

1 -0.8 


4 


9+- 


3.74c-14 


30 


37 


HLP-1,B89-2,0541458-0154297,B-92,M-205,V- 


188*f 


J054145. 8-015410 


5 


41 


45.89 


-1 


54 


10.9 


1065 ± 34 


2.53 


3.31 


3 


n+0.4 
y -0.3 


2 


0+0:5 

— U.4 


2.47e-13 


31 


25 


M-206,V-20 S 


189* 


J054145.9-015626 


5 


41 


45.94 


-1 


56 


26.4 


176 ± 14 


1.45 


3.88 


5 


c+2.6 

°-1.9 






5.03c-14 


30 


54 


05414594-0156263,M-207,B89t23 


190 


J054145. 9-015501 


5 


41 


45.97 


-1 


55 


01.9 


65 ± 9 


1.61 


3.41 










2.42c-14 






HLL-96,HLP-92,M-208 


191 


J054146. 0-015653 


5 


41 


46.02 


-1 


56 


53.8 


33 ± 7 


1.38 


3.76 










8.40e-15 






05414602-0156538 


192*f 


J054146. 1-015414 


5 


41 


46.11 


-1 


54 


14.7 


1753 ± 43 


27.4 


3.89 


3 


7+0.5 

'-0.5 






5.20e-13 


var. 


05414611-0154147,B-93,M-210 


193* 


J054146.1-015621 


5 


41 


46.15 


-1 


56 


21.9 


397 ± 21 


2.96 


4.24 


9 


6 +2.4 
D -1.8 


5 


o+H 

— i.y 


1.27c-13 


31 


07 


V-21 


194 


J054146.2-015346 


5 


41 


46.20 


-1 


53 


46.7 


66 ± 9 


0.79 


2.49 






2.00c-14 






HLL-58,05414621-0153466,M-212,B89-7 


195 


J054146.2-015533 


5 


41 


46.23 


-1 


55 


33.1 


165 ± 14 


1.89 


4.15 


11.811°;° 


3 


o+... 

°-1.9 


4.84e-14 


30 


86 


M-211 


196* 


J054146.2-015654 


5 


41 


46.23 


-1 


56 


54.6 


150 ± 13 


1.17 


4.07 


5 


9+2.4 






4.50e-14 


30 


41 


V-22 


197 


J054146. 2-015554 


5 


41 


46.26 


-1 


55 


54.1 


32 ± 7 


0.34 


4.77 










1.18c-14 








198 


J054146. 4-015011 


5 


41 


46.40 


-1 


50 


11.3 


47 ± 8 


0.91 


2.01 










5.64e-15 






HLL-5,05414641-0150111 


199f 


J054146. 5-015148 


5 


41 


46.52 


-1 


51 


48.6 


91 ± 11 


6.37 


3.42 


3 


n +1.3 
U -0.5 


4 




2.33c-14 


30 


17 


05414652-0151487 


200 


J054146.5-015446 


5 


41 


46.57 


-1 


54 


46.6 


94 ± 11 


1.24 


3.27 


7 


n +5.7 
u -2.2 


1 




6.68e-14 


31 


06 


05414655-0154469,B-94,M-214,V-23 


201 


J054146. 8-015447 


5 


41 


46.86 


-1 


54 


47.4 


16 ± 5 


1.26 


4.14 










2.70c-14 






M-217 


202 


J054146. 8-014957 


5 


41 


46.89 


-1 


49 


57.1 


370 ± 20 


0.88 


2.07 





O+0.2 
S -0.2 


2 


c+0.7 

°-0.6 


4.48e-14 


30 


25 


HLL-4,HLP-23,05414690-0149573 


203 


J054147.0-014926 


5 


41 


47.01 


-1 


49 


26.6 


449 ± 22 


3.14 


2.47 


1 


c+0.3 
°-0.4 


2 


9+1.5 
y -0.7 


7.13c-14 


30 


55 


HLL-1, 05414702-0149265 


204 


J054147.0-015035 


5 


41 


47.10 


-1 


50 


35.4 


158 ± 14 


0.51 


1.93 


1 


+0.4 
u -0.4 


1 


O + 0.4 

•^-o.s 


1.67c-14 


29 


99 


HLL-8,HLP-90,05414712-0150353 


205 


J054147.1-020018 


5 


41 


47.14 


-2 


00 


18.8 


29 ± 7 


0.96 


1.65 








2.27c-15 






05414714-0200190 


206 


J054147.2-015819 


5 


41 


47.21 


-1 


58 


19.7 


12 ± 5 


0.94 


2.81 










1.98c-15 






HLL-138,05414723-0158202 


207*f 


J054147.4-015526 


5 


41 


47.45 


-1 


55 


26.3 


899 ±31 


17.9 


3.74 


4 


,+0.7 
°-0.6 


7 


8 -2.8 


2.40c-13 


31 


11 


M-220 


208 


J054147.9-015541 


5 


41 


47.93 


-1 


55 


41.3 


24 ± 6 


0.55 


4.06 










6.80c-15 






05414793-0155416,M-225 


209 


J054147.9-015358 


5 


41 


47.97 


-1 


53 


58.4 


17 ± 5 


0.67 


3.45 










4.13c-15 






05414796-0153585.M-226 
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No. 


Name 


R.A. 


Dccl. 


Net 


KS 


E 


N H 


kT 


Fx,abs 


log L x 


Identification 






(J2000) 


(J2000) 


Counts 




(keV) 


(10 22 cm" 2 ) 


(kcV) 


(erg/cm 2 /s) 


(erg/s) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 



210*f 


J054148. 2-015601 


5 


41 


48.21 


-1 


56 


01.8 


1616 ± 41 


15.4 


3.61 


4.5+°-? 

—0.3 


—0.9 


4.17c-13 


31.45 


05414821-0156020,M-230,V-24 


211 


J054148. 6-015349 


5 


41 


48.60 


-1 


53 


49.5 


95 ± 11 


0.96 


2.52 


—0.9 


3.0+V q 

— 1.3 


1.50e-14 


29.80 


HLL-59,05414861-0153495 


212 


J054148. 6-015540 


5 


11 


48.63 


-1 


55 


41.0 


9 ± 4 


0.99 


4.65 






3.49e-15 






213 


J054148. 6-015415 


5 


41 


48.69 


-1 


54 


15.6 


32 ± 7 


2.42 


2.59 






5.46c-15 




HLL- 74, 05414868-0154159 


214 


J054149. 3-015332 


5 


41 


49.31 


-1 


53 


32.7 


76 ± 10 


1.94 


2.66 


1.0 


3 l + ;' v. 
Ji — 1.6 


1.25e-14 


29.79 


HLL-50, 05414934-0153318 


215 


J054149. 6-014601 


5 


41 


49.63 


_1 


46 


01.7 


146 ± 13 


2.34 


2.03 


x j+0.5 
1- -0.4 


1 6 +0 ' 7 


1.70e-14 


30.05 


05414961-0146015 


216 




5 


41 


49.66 


_l 


53 


26.9 


89 i 10 


0.99 


2.79 


1 6+ ' 7 
1 '°-0.6 


4 8~*~"" 
4 -°-2.3 


1.53e-14 


29.78 


HTJ, 48 r >414Qfifi D1 ^971 


217 


T0541 49 8-01 5255 


5 


41 


49.81 


_l 


52 


55.8 


21 ± 6 


0.98 


4.37 




6.63e-15 




HT,T,-36 05414983-0152560 


218 


T0541 4Q 8-01 5947 


5 


41 


49.88 


_l 


59 


47.4 


11 ± 5 


0.47 


3.39 






2 76e-1 5 




0541 4984-01 59475 


219 


T054150 0-015807 


5 


41 


50.04 


-1 


58 


07.8 


11 ± 4 


1.29 


1.56 






9.11e-16 




0541 5006-01 58083 


220 


J054150. 1-015744 


5 


41 


50.15 


-1 


57 


44.7 


154 ± 13 


0.95 


3.12 


2 5+ 1 ' 


4 2+ 1 ' 8 
-1.8 


3.26c-14 


30.18 


HLL-135 HLP-98 05415014-0157449 


221 


T0541 50 5-01 531 5 


5 


41 


50.50 


_x 


53 


15.2 


kk8 + 25 

•J 'JO _l_ iiij 


1.25 


2.11 


9 +0 ' 2 

u - a -o.i 


9 n+0.4 
z - -0.3 


6.78e-14 


30.48 


HT,T,-41 0541 5050-01 531 53 R89-6 


222*f 


10^41^0 ^ 0909^fi 


5 


11 


50.50 


_2 


02 


56.9 


242 ± 17 


6.10 


3.46 


4 5+ 1 ' 3 
4 -°-0.9 


9+2.5 


6.34e-14 


30.67 


0^41^0^^ 0909^70 

UiJ4:lUUJO -UZUZO I u 


223* 


to^41 p. 01 ^nfi 


r: 



41 




_1 






99 -4- 
ZZ HI u 


n R8 
u.uo 


1 8fi 






9 n^fj 1 ^ 




sec notGS — 


224 


T0^41 ^fl 7 01 ^48 


5 


11 


50.74 


_l 


55 


48.3 


13 i 5 


0.75 


2.76 












225 


T0^41 ^fl 8 01 ^1 9 


5 


41 


50.86 


_l 




12.7 


22 i 6 


0.77 


2.67 






^ 40p 1 ^ 




HTJ, 40 0^41 ^0S7 01 ^^1 9K 


226 


10^41 & m ^799 


r: 



41 


87 


_l 


n.7 


22 2 






a 






9 fi^o 1 ^ 






227 


10^41 ^0 q m 4&01 


r: 



11 


ou.yo 


_l 


18 


01 




n 7 a 


o.uo 






9 ^Qr> 1 ^ 






228 


to^41 ^0 q m ^90 


r: 



11 


^1 on 


_l 




21 


1 ^ -+- ^ 


u.oy 


1 1(1 






94^ 1 fi 
. z 10- 1 u 




0^41 ^1 00 01 ^^91 ^ 


229 


J054151. 6-015216 


5 


41 


51.60 


-1 


52 


16.6 


66 ± 9 


0.50 


2.34 






1.00c- 14 




05415161-0152167 


230 


J054152. 0-015518 


5 


41 


52.08 


-1 


55 


18.1 


18 ± 5 


1.16 


1.58 






1.58e-15 




05415209-0155183 


231 


J054152. 1-015717 


5 


41 


52.15 


-1 


57 


17.3 


30 ± 7 


1.37 


2.34 






1.35c-14 




HLL-127,HLP-75,B89-26,05415215-0157174 


232* 


J054152. 1-020121 


5 


41 


52.20 


-2 


01 


21.1 


171 ± 14 


0.77 


2.21 


1 r+0.6 
1 '°-0.B 


1 6+°' 6 


2.41e-14 


30.17 


05415220-0201211 


233 


J054152.9-015629 


5 


41 


52.93 


-1 


56 


29.3 


55 ± 8 


0.53 


3.35 






1.42e-14 




HLL-121,B89-25,05415292-0156297 


234*f 


J054152. 9-015634 


5 


41 


52.94 


-1 


56 


34.8 


783 ± 29 


6.04 


2.58 


1 r+0.3 
^-O^ 


3 1+ 11 


1.31c-13 


30.80 


HLL-123,05415294-0156351 


235 


J054153. 1-015743 


5 


41 


53.19 


-1 


57 


43.8 


59 ± 9 


1.91 


3.57 






2.96e-14 




HLL-134,HLP-78,05415319-0157441,B89-27 


236 


J054153.8-015516 


5 


41 


53.87 


-1 


55 


16.1 


304 ± 18 


1.04 


2.77 


2 5+ ' 7 


9 n+0.8 


5.19c-14 


30.55 


HLL-101,HLP-88,05415387-0155163,B89-21 


237 


J054154. 1-015817 


5 


41 


54.12 


-1 


58 


17.1 


12 ± 5 


0.59 


1.50 






1.57c-15 




05415406-0158170 


238 


J054154.1-014801 


5 


41 


54.19 


-1 


48 


01.7 


77 ± 10 


1.26 


1.93 






1.02e-14 




05415419-0148017 


239 


J054154.3-015556 


5 


41 


54.31 


-1 


55 


56.2 


25 ± 6 


1.01 


1.22 






1.67c-15 




HLL-112,05415430-0155566 


240 


J054154.7-015726 


5 


41 


54.70 


-1 


57 


27.0 


33 ± 7 


0.49 


1.61 






2.74c-15 




HLL-128,05415472-0157273 


241 


J054154.9-015647 


5 


41 


54.93 


-1 


56 


47.9 


22 ± 6 


0.56 


2.24 






3.16e-15 




HLL-125,05415492-0156483 


242 


J054154.9-015630 


5 


41 


54.96 


-1 


56 


30.9 


18 ± 5 


1.20 


3.22 






3.82e-15 




05415497-0156312 


243 


J054154.9-015957 


5 


41 


54.99 


-1 


59 


57.8 


19 ± 6 


0.67 


3.72 






5.00c-15 






244 


J054155. 0-020030 


5 


41 


55.04 


-2 


00 


30.9 


111 ± 12 


3.36 


3.13 


4 2+ 1 ' 6 


2 5+ 1 ' 4 
z -°-0.9 


2.34c-14 


30.30 


05415505-0200309 



Table 1 — Continued 



No. 


Name 


R.A. 


Decl. 


Net 


KS 


E 


Nh 


kT 


1 x,abs 


log Lx 


Identification 






(J2000) 


(J2000) 


Counts 




(kcVI 


(10 22 cm- 2 ) 


(kcVI 


(OTP / CTTI^ /s^ 


('ore' /si 




(1) 


(21 


(31 




(4) 


(51 


(61 


(71 


(81 


(91 


(101 


fill 
V- 1 - 


(12) 


245 


J054155.2-014847 


5 41 55.21 


-1 


48 47.4 


15 ± 6 


0.75 


1.85 






1.21e-15 




05415518-0148485 


246 


J054155.7-015508 


5 41 55.79 


-1 


55 08.5 


135 ± 13 


1.61 


1.40 


o"i±g-? 


1 4+0-5 

— 0.7 


9.95c-15 


29.41 


HLL- 99,05415579-0155092 


247 


J054156.5-014616 


5 41 56.57 


-1 


46 17.0 


217 ± 16 


0.82 


2.11 


6 +0 ' 3 

U.D_g 2 


2 ^ +1 - 5 

7 


2.78c-14 


30.17 


05415656-0146168 


248 


J054156.6-015552 


5 41 56.68 


-1 


55 52.4 


170 ± 14 


3.12 


1.78 


0.b+_° % 


i.5±°;« 


1.56e-14 


29.78 


HLL-110,05415668-0155526 


249 


J054156.7-015853 


5 41 56.74 


-1 


58 53.4 


28 ± 6 


0.49 


2.29 


4.33c-15 




05415674-0158538 


250 


J054156.7-015351 


5 41 56.79 


-1 


53 51.8 


264 ± 17 


1.92 


2.77 


1.7 +0 l 

— U.4 


3.6± 2 ;2 


1.18c-13 


30.71 


05415680-0153521 


251 


J054157.0-015032 


5 41 57.02 


-1 


50 32.0 


16 ± 5 


1.47 


1.86 






1.61e-15 




HLL-7,05415697-0150331 


252 


J054157.8-015127 


5 41 57.83 


-1 


51 27.7 


38 ± 7 


1.22 


1.49 






2.64e-15 




HLL-14,05415786-0151278 


253 


J054158.0-015957 


5 41 58.10 


-1 


59 57.3 


39 ± 7 


1.82 


1.37 






3.09c-15 




05415807-0159575 


254 


J054159.8-014735 


5 41 59.82 


-1 


47 35.3 


39 ± 8 


1.29 


2.04 






3.23c-15 




05415976-0147363 


255 


J054159.8-015949 


5 41 59.89 


-1 


59 49.6 


7 ± 4 


0.79 


3.49 






2.20c-15 




05415997-0159496 


256 


J054201. 1-015621 


5 42 01.14 


-1 


56 21.3 


23 ± 6 


0.63 


1.70 






2.66e-15 




05420114-0156214 


257 


J054201.4-020008 


5 42 01.47 


-2 


00 08.3 


15 ± 5 


0.36 


3.74 






3.77e-15 






258* 


J054201. 7-015629 


5 42 01.71 


-1 


56 29.6 


160 ± 14 


1.84 


2.29 


9+ - 3 


o 1+2.4 
— 1 1 


2.19c-14 


29.94 


HLL-120,05420172-0156299 


259 


J054201. 9-015629 


5 42 01.98 


-1 


56 29.2 


19 ± 5 


0.48 


1.77 




1.64e-15 




05420201-0156289 


260 


J054204.2-015201 


5 42 04.27 


-1 


52 01.8 


16 ± 5 


0.66 


3.21 






3.40e-15 






261 


J054205.2-015344 


5 42 05.27 


-1 


53 44.4 


65 ± 9 


2.38 


1.87 






9.36e-15 




05420528-0153447 


262 


J054205.7-015928 


5 42 05.77 


-1 


59 28.7 


42 ± 8 


1.46 


2.49 






7.67c-15 




05420580-0159289 


263 


J054205.9-015252 


5 42 05.95 


-1 


52 52.0 


136 ± 13 


1.53 


1.79 


4+ ' 2 


n 9+4.1 
z ' z -0.6 


3.63c-14 


30.08 


HLL-35,05420597-0152522 


264 


J054207.5-015335 


5 42 07.51 


-1 


53 35.1 


22 ± 6 


0.52 


1.75 






1.68G-15 




05420744-0153368 


265 


J054207.5-015738 


5 42 07.55 


-1 


57 38.2 


18 ± 6 


0.72 


3.66 






6.05e-15 






266 


J054207.6-015322 


5 42 07.60 


-1 


53 22.1 


9 ± 4 


0.52 


3.05 






1.73c-15 






267 


J054208.9-015741 


5 42 08.98 


-1 


57 41.2 


25 ± 6 


0.77 


2.83 






4.68c-15 






268 


J054209.3-015804 


5 42 09.35 


-1 


58 04.5 


528 ± 24 


1.37 


1.78 


0.5+°- 2 


2.0+0-1 


5.32c-14 


30.27 


05420934-0158048 


269* 


J054209.5-020004 


5 42 09.57 


-2 


00 04.2 


1207 ± 36 


4.47 


1.95 




9 4 +0.3 


1.36e-13 


30.69 


05420956-0200044 


270 


J054210.2-020116 


5 42 10.22 


-2 


01 16.3 


10 ± 5 


0.72 


2.74 






1.25e-15 






271 


J054211.0-014626 


5 42 11.07 


-1 


46 26.2 


24 ± 6 


0.76 


2.09 






2.99e-15 




05421113-0146255 


272* 


J054211. 2-015943 


5 42 11.25 


-1 


59 43.3 


1352 ± 38 


1.53 


1.26 


0.2t0; 2 




1.18e-13 




05421125-0159437 


273 


J054211. 8-015637 


5 42 11.89 


-1 


56 37.9 


9 ± 4 


0.98 


2.86 






1.78c-15 






274* 


J054213.1-015833 


5 42 13.19 


-1 


58 33.3 


1251 ± 36 


2.85 


2.07 


7+°- 2 
u -'-0.2 


9 c+0.4 


1.72o-13 


30.81 


05421319-0158336 


275 


J054214.1-015141 


5 42 14.12 


-1 


51 41.3 


13 ± 6 


0.53 


4.17 


6.90o-15 






276 


J054215.5-015500 


5 42 15.55 


-1 


55 00.1 


15 ± 5 


0.96 


3.07 






3. 03c- 15 






277 


J054216.2-015614 


5 42 16.24 


-1 


56 14.6 


180 ± 15 


1.45 


1.40 


2+°- 2 


o-8±8:i 


1.44e-14 


29.80 


Haro 5-53,05421620-0156148 


278 


J054216.7-015308 


5 42 16.79 


-1 


53 08.9 


19 ± 6 


0.99 


2.10 






2.80c-15 




05421677-0153091 


279 


J054218.0-015127 


5 42 18.08 


-1 


51 27.9 


20 ± 6 


0.77 


2.78 






4.74e-15 







GO 



Table 1 — Continued 



No. 
(1) 


Name 
(2) 


R.A. 
(J2000) 
(3) 


Dccl. 
(J2000) 

(4) 


Net 
Counts 
(5) 


KS 

(6) 


E 

(keV) 
(7) 


N H 
(10 22 cm~ 2 ) 
(8) 


kT 

(keV) 
(9) 


Fx.abs 

(erg/cm 2 /s) 
(10) 


log L x 
(erg/s) 
(11) 


Identification 
(12) 


280 
281 
282 
283 


J054218.4-015731 
J054218.6-015421 
J054221. 2-015910 
J054226.5-015913 


5 42 18.45 
5 42 18.63 
5 42 21.29 
5 42 26.53 


-1 57 31.7 
-1 54 21.8 
-1 59 10.2 
-1 59 13.2 


77 ± 10 
16 ± 6 
88 ± 11 
8 ± 5 


0.72 
0.60 
0.82 
1.44 


3.39 
3.24 
1.64 
2.63 


O.ot ; 1 


1 0+ ' 3 
1 - u -0.2 


2.31e-14 
5.04e-15 
l.Olc-14 
2.19e-15 


29.75 


05422123-0159104 



a Col. (1) running source number; an asterisk indicates that additional source notes follow below; an / inidicatcs that a flare or other obvious variability is 
seen in the source light curve, (2) J2000 IAU source name, (3)-(4) J2000 Chandra position with astrometric correction applied, (5) net source counts (= total 
counts - background counts) ± net counts uncertainty (Gchrels 1986) in the 0.5 - 7.0 keV range; source counts were measured in PSF-corrected 95% encircled 
energy regions except for the brightest sources where 99% encircled energy regions were used; smaller regions were used for some closely-spaced sources, (6) 
KS variability statistic, sources with KS > 1.7 are variable (see text), (7) mean photon energy as determined from source events in the 0.5 - 7 keV range, 
(8)-(9) absorption column density N// and time-averaged characteristic temperature kT for sources with >90 counts as determined from spectral fits with IT 
VAPEC optically thin plasma model using a fixed iron abundance Fe = 0.3 solar (errors are 90% confidence), (10) absorbed flux (0.5 - 7 keV) as determined 
from source event list (see text), (11) unabsorbed luminosity (0.5 - 7 keV) at d = 415 pc determined from spectral fits of brighter sources where an estimate 
of N H and kT is available, (12) candidate IR or optical identification from 2MASS, SIMBAD, and B89 = Barnes et al. 1989; B = Beck et al. 2003; HLL = 
Haisch, Lada & Lada 2000; HLP = Haisch et al. 2001, M = Meyer 1996. IR identifications found in the all-sky release of the 2MASS Point Source Catalog 
are listed by their coordinates (e.g. 05411324-0153306). Radio identifications (V-xx) are from 3.6 cm VLA observations of Rodriguez et al. 2003. 
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5. IT model gives a high but uncertain temperature kT > 4 keV. 

15. HLL-131 is a likely IR counterpart, but there is a ~4" positional offset even after correcting for system- 
atic differences between CXO and HLL coordinates. 

16. HLP-54 is a likely IR counterpart, but there is a w4" positional offset between CXO and HLP coordi- 
nates. 

18. Flare spectrum gives kT = 5.4+ 2 " 2 keV and log Lx = 31.58 ergs s^ 1 (unabsorbed). 

27. Flare spectrum gives kT = 4.2^ keV and log Lx = 30.99 ergs s" 1 (unabsorbed). (Fig. 7) 

36. Variability is present throughout the observation. Complex or rapidly varying temperature structure 

suspected. 

39. Source lies in CCD gap. Flux is uncertain. 

41. Weak flare whose spectrum gives kT = 2.6^0.5 keV and log Lx = 30.80 ergs s _1 (unabsorbed). 
52. Weak flare. Insufficient counts to accurately determine flare parameters. 

75. IT models give Xred = 1-09 but 2T models give a slightly better fit with Xred = 1-05 and converge to Nh 
= l-7±g"J x 10 22 cm" 2 , kTx = 0.6t Q i keV, kT 2 = 2.7±° j keV, and log Lx = 31.50 ergs s" 1 (unabsorbed). 
81. Flare spectrum gives kT = 9.8t 4 " 5 keV and log Lx = 31.08 ergs s _1 (unabsorbed). (Fig. 7) 
85. Flare spectrum gives kT = 3.8±g;| keV and log Lx = 31.34 ergs s" 1 (unabsorbed). (Figs. 7 & 10) 

93. Flare spectrum gives kT = 3.5± 2 ;i keV and log Lx = 30.98 ergs s (unabsorbed). 

94. Slow variability is present throughout the observation. DEM model shows a double-peaked structure 
with cool and hot components. 

99. IT models give xfed ~ 0-94 but 2T models give a slightly better fit with xfed = 0-92 and converge to a 
larger absorption N H = 2.7t°;g x 10 22 cm -2 . 

104. Declining count rate may indicate decay phase of a flare. 

105. Weak flare. Insufficient counts to accurately determine flare parameters. 

108. Flare spectrum gives same temperature kT = 3.3^ 'g keV and log Lx = 31.01 ergs s -1 (unabsorbed). 
114. Declining count rate may indicate decay phase of a flare. Possible spectral contamination from nearby 
source. 

120. Source may be affected by moderate pilcup (wlO - 15%) during short periods when the count rate 
fluctuated to levels of «0.15 c s~ 1 . 

123. Count rate increased throughout the observation (Fig. 7). 

133. IT models are only marginally acceptable with xfed = 1-35. DEM model peaks near 1.7 keV but shows 

substantial emission up to ~4 keV. Declining count rate may indicate decay phase of a flare. 

139. Slowly rising variability and outburst spectrum gives kT = 4.2^'g keV and log Lx = 30.79 ergs s -1 

(unabsorbed). 

141. The spectrum used to determine Nh, kT excludes data from a weak flare that occurred during the last 
5 ksec of the observation. 

145. No reliable spectral fit because of contamination from nearby source no. 148. Possible optical counter- 
part is GSC 04771-00589. 

148. No reliable spectral fit because of contamination from nearby source no. 145. 

152. IT models are only marginally acceptable with xfed = 1-38. DEM model peaks near 1.4 keV. 

153. IT models are not acceptable with xf e d = 1-91 but 2T models are acceptable with xf e d = 1.00. Nr is 
from 2T model (kTi = 0.8 keV, kT 2 = 3.4 keV). 

160. Weak flare. Insufficient counts to accurately determine flare parameters. 

174. Flare spectrum gives kT = 4.5^7 keV and log Lx = 31.45 ergs s -1 (unabsorbed). Moderate pilcup 
may have occurred near flare peak. (Fig. 7) 

175. Long duration flare with a secondary peak during decay phase (Fig. 7). Flare spectrum gives kT = 
5.2t 2 ;g keV and log Lx = 31.47 ergs s" 1 (unabsorbed). 



-61 - 



182. The counterpart name IRS 2b is from Nisini et al. (1994). Spectral fit parameters are uncertain because 
of possible contamination from source 183 located 1.7" to the north (Fig. 2). Some improvement is obtained 
with 2T models which place most of the emission measure in a cool component at kTi ss 0.2 keV and require 
a larger absorption column N# = 5.5+ 2 'q x 10 22 cm ~ 2 - 

183. Weak flare. Insufficient counts to accurately determine flare parameters. Spectral fit parameters are 
uncertain due to possible contamination from source 182. 

187. IT models give xl e d = 1-10 ^T models give a slightly better fit with xf e d = 0-88 and converge to a 
higher absorption N# = 6.1^7 x 10 22 cm~ 2 , kTi = 0.6jz 3 keV, and a hotter component with uncertain 
temperature kT2 > 3 keV. 

188. Weak flare. Flare spectrum gives kT = 3.2+{j;g keV and same log Lx = 31.25 ergs s" 1 (unabsorbcd). 

189. IT models are not acceptable (x^ed = 1-52) nor are 2T models (xled = 1-68). Hard heavily-absorbed 
spectrum with uncertain temperature. 

192. Slowly rising outburst (Fig. 7) with high absorption. The X-ray temperature shortly after onset reached 
values of at least kT ~ 6 keV, and possibly as high as kT ~ 10 keV. The unabsorbed luminosity varied from 
a pre-outburst value log Lx = 30.5 ergs s _1 to at least log Lx = 31.8 ergs s _1 during the outburst. Source 
is affected by moderate pileup (w27%) near end of observation. 

193. Declining light curve may indicate decay phase of a flare. 

196. IT models give a high but uncertain temperature kT > 3.1 keV. 

207. Slowly rising outburst (Fig. 7). Very hard source. Time-averaged outburst spectrum gives kT = 7.8+3.3 
keV and log Lx = 31.4 ergs s _1 (unabsorbed). Temperature peaked about 4 hours after onset followed by 
a general decline, but there is evidence for at least one reheating event. DEM model of outburst spectrum 
suggests that most of the emission measure comes from plasma above ~4 keV. Insufficient counts to analyze 
pre-outburst spectrum. 

210. Declining light curve during first ^20 ksec of observation may indicate decay phase of a flare. Decay 
spectrum gives kT = 4.9+J'I keV and log Lx = 31.59 ergs s _1 (unabsorbed). Source may be affected by 
moderate pileup («16%) at beginning of observation. 

222. Weak flare. Insufficient counts to accurately determine flare parameters. 

223. A faint IR source is visible at the X-ray position in 2MASS H and K-band images (2MASS J05415068- 
0153063), but is flagged as a filter glint artifact. Because of the close positional agreement with Chandra, 
the IR source flagged as an artifact may be real. 

232. IT models are only marginally acceptable with Xred = 1-39. 
234. Declining light curve may indicate decay phase of a flare. 
258. Possible double source. 

269. Count rate declined throughout the observation (Fig. 7). 

272. Soft source (low absorption). IT models underestimate the flux below 1 keV and are not acceptable 
(Xred = 3.2). 2T models give improved fits when the Fe abundance is allowed to vary (x 2 ed = 1-24), con- 
verging to Fe = 0.8 solar with kTi = 0.77 keV and kT 2 = 1.9 keV. DEM model peaks near 0.7 keV. 
274. IT models are only marginally acceptable with x 2 re d = 1-42- DEM model peaks near 2 keV. Variability 
is present throughout the observation. 



